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SECTION  1 


INTRODUCTION 


1 .1  Scope 

Volume  2  of  the  present  report  giver  an  account  cf  theoretical 
advances  and  associated  demonstrations  made  during  the  reporting  period 
relating  to  the  synthesis  of  active  controllers  for  large  flexible 
structures.  Particular  applications  discussed  include  the  quenching  of 
vibrations  in  support  structures  for  optical  pointing  systems  and  the 
efficient  generation  of  closed-loop  large-angle  slewing  maneuvers  for 
linear  tracking.  Certain  aspects  of  discrete-time  dynamic  modeling  for 
structural  vibration  control  are  also  treated  in  depth.  The  vehicles  for 
the  demonstration  examples  are  the  ACOSS  Model  No.  2  representation  of  an 
optical  pointing  system  (for  vibration  control)  and  a  rotating  rigid  hub 
with  four  flexible  appendages  (for  large-angle  slew  control).  A  concise 
summary  of  the  principal  results  contained  in  the  present  volume  is  given 
in  the  remainder  of  Section  1 . 

1.1.1  Vibration  Control 

Effective  accommodation  of  disturbances  is  fundamental  to  closed- 
loop  regulation,  being  especially  vital  for  the  quenching  of  vibrations 
in  large  flexible  optical  pointing  systems.  Accommodation  of  a  (small) 
finite  number  of  periodic  disturbances  at  fixed  known  frequencies  poses 
only  mild  difficulty.  Substantially  more  difficult  is  the  problem  of 
accommodating  an  aperiodic  disturbance  having  nonnegligifcle  power 
spectral  density  across  a  broad  range  of  frequencies.  In  the  preceding 
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reporting  period,  a  scientific  experiment  was  conceived  to  examine  the 
question  of  active  control  synthesis  for  broadband  disturbance  accommoda¬ 
tion  in  which  the  selection  of  reduced-order  models  and  of  active 
transducers  were  to  be  given  substantial  consideration. 

Sections  2  and  3  of  the  present  volume  taken  together  constitute 
an  initial  account  of  progress  in  the  study  of  broadband  disturbance 
accommodation  for  large  flexible  optical  pointing  systems.  The 
disturbance  to  be  neutralized  has  a  constant  nonzero  power  spectral 
density  over  the  frequency  range  between  0  and  5  Hz.  In  the  first  part 
(Section  2),  attention  is  focused  primarily  upon  the  selection  of  active 
transducers,  both  by  functional  type  and  by  location.  The  theoretical 
basis  of  a  recently-developed  selection  algorithm  is  fully  described;  the 
principal  elements  are  full-rank  matrix  factorizations  --  in  particular, 
the  QR-decomposition  —  and  their  employment  in  the  representation  of 
solutions  to  least-squares  approximation  problems.  The  algorithm  is 
designed  to  select  a  relatively  small  number  of  actuators  that  have  a 
relatively  high  influence  on  the  variables  to  be  regulated  (in  this  case, 
optical  system  pointing  errors).  Results  of  four  examples  with  ACOSS 
-fodel  No.  2  demonstrate  that  selections  are  readily  achievable  that  are 
substantially  fewer  in  number  than  the  number  of  controlled  modes,  and 
yet  give  promise  of  effective  performance.  In  the  second  part  (Section 
3),  attention  is  shifted  to  the  design  and  evaluation  of  controllers  for 
ACOSS  Model  No.  2  associated  with  the  four  active  transducer  selections 
described  in  Section  2.  The  controller  feedback  structure  is  generated 
using  a  textbook  disturbance-rejection  control  design.  Attempts  at 
stability-enhancing  adjustments  beyond  meeting  the  specifications  on  the 
optical  pointing  error  are  deliberately  postponed  in  order  to  isolate  the 
influence  upon  the  overall  synthesis- process  of  the  selection  of  reduced- 
order  models  and  of  active  transducers,  including  their  mutual  inter¬ 
action.  In  each  example,  the  effect  of  unmodeled  modes  upon  control 
system  scability  and  performance  is  examined  using  an  expanding  family  of 
reduced-order  evaluation  models.  In  the  evaluations,  closed-loop 
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stability  over  a  range  of  frequencies  substantially  beyond  that  of  the 
design  model  is  demonstrated. 

In  the  modeling  of  dynamic  systems,  state  space  representation 
provides  a  unified  framework  which  abstracts  certain  system-theoretic 
characteristics  common  to  a  wide  variety  of  physical  systems.  The 
resulting  separation  of  dynamic  characteristics  from  specific  physical 
processes  provides  beneficial  clarification.  However,  it  has  the 
disadvantage  of  masking  the  physics  of  the  problem.  The  sharpened 
understanding  that  follows  from  taking  into  account  the  unique 
characteristics  of  a  specific  physical  process  usually  increases  the 
performance  achievable  from  a  process  controller.  In  this  spirit,  a 
frrsh  approach  to  sampled-data  control  for  flexible  structures  is  taken 
in  Section  4  i-  which  second-order  difference  equations  provide  the 
dynamic  model.  For  the  problem  of  velocity  output  feedback,  explicit 
sufficient  conditions  for  stability  in  terms  of  the  feedback  gain  matrix 
and  the  sampling  rate  are  given.  Surprisingly,  these  conditions  require 
non-colocation  of  actuators  and  sensors  as  a  function  of  the  sampling 
rate.  The  non-colocation  property  collapses  to  the  familiar  colocation 
property  for  analog  systems  as  the  sampling  period  converges  to  zero. 

1.1.2  Large-Angle  Slew  Control 

Large-angle  slewing  maneuvers  for  flexible  spacecraft  generally 
require  the  consideration  of  nonlinear  structural  dynamics.  However 
certain  problems  of  substantial  interest  for  applications  require  only 
the  consideration  of  linear  autonomous  structural  dynamics.  Several  such 
problems  are  formulated  and  solved  in  Section  5;  namely,  the  problem  of 
optimally  tracking  a  known  trajectory  (without  hard  terminal 
constraints),  and  the  problem  of  optimal  linear  slew  in  which  the 
terminal  values  of  some  or  all  of  the  state  variables  at  a  fixed  finite 
terminal  time  are  prescribed.  The  emphasis  is  upon  obtaining  explicit 
expressions  for  solutions  of  the  necessary  conditions  for  optimality  in 
terms  of  matrix  parameters  which  are  solutions  of  linear  algebraic 
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equations.  The  problem  of  computing  solutions  is  thereby  reduced  from 
one  of  numerical  calculus  (i.e. ,  integration)  to  one  of  numerical  linear 
algebra.  Significant  increases  in  computational  efficiency  resulting 
from  this  approach  are  demonstrated  with  a  variety  of  examples  on  a 
rotating  rigid  hub  with  four  flexible  appendages.  The  method  of 
obtaining  explicit  solution  expressions  is  also  shown  to  have  direct 
application  to  the  efficient  generation  of  feedback  control  solutions  to 
the  general  unconstrained  fixed  time  linear  optimal  control  problem. 

1 .2  Limitations 

The  scientific  experiment  relating  to  active  control  synthesis  for 
broadband  disturbance  ac<"  .nmodation  is  incomplete.  Projected  work 
includes: 

(1)  Expansion  of  the  class  of  candidates  from  which  active 
transducers  are  selected; 

(2)  Examination  of  the  interaction  between  the  selection  of 
reduced-order  models  and  of  active  transducers;  and 

(3)  Modifications  of  the  controller  feedback  structure  to  ensure 
an  acceptable  compromise  between  the  requirement  to 
neutralize  the  effects  of  disturbances,  and  the  limitations 
imposed  by  an  incomplete  knowledge  of  structural  and 
disturbance  characteristics. 
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PROGRESS  ON  SYNTHESIS  OF  ACTIVE  CONTROL  FOR  BROADBAND 
DISTURBANCE  ACCOMMODATION 

PART  1:  ACTIVE  TRANSDUCER  SELECTION 


2.1  Background 

In  the  automatic  control  of  dynamic  systems,  the  achievement  of 
stability  and  a  specified  lev«i  of  performance,  in  spite  of  anticipated 
disturbances,  is  fundamental.  When  the  object  of  control  has  substantial 
inherent  flexibility,  the  problem  of  accommodating  disturbances  is 
aggravated  by  the  tendency  of  the  disturbances  to  excite  structural 
vibrations.  Even  with  flexible  structures,  the  effects  of .a  few  periodic 
disturbances  at  known  frequencies  are  relatively  easy  to  accommodate. 
However,  when  the  frequency  of  periodic  disturbances  is  unknown  or 
variable,  or  when  the  disturbances  are  aperiodic,  disturbance 
accommodation  i3  much  more  difficult.  Relatively  little  has  been  written 
regarding  accommodation  of  the  latter  class  of  disturbances  in  the 
context  of  large  flexible  structures.  One  careful  study  treating 
unknown-frequency  periodic  disturbances  [Ref.  2-1}  has  reported  moderate 
success.  It  is  notable  that  the  treatment  reported  gave  careful 
attention  to  the  placement  of  active  transducers  (i.e. ,  actuators  and 
sensors).  The  only  study  available  to  date  on  accommodation  of  aperiodic 
disturbances  with  a  continuum  frequency  spectrum  ("broadband")  in  the 
context  of  large  flexible  structures  [Ref.  2-2)  also  reports  moderate 
success,  without  however,  giving  any  attention  to  the  selection  of  active 


transducers  beyond  the  analysis  developed  for  accommodating  periodic 
disturbances  of  known  frequencies  [Ref.  2-3] . 

The  presence  of  a  broadband  disturbance  escalates  the  importance 
of  giving  careful  attention  to  the  selection  of.  reduced-order  models,  the 
selection  of  active  transducers,  and  the  interdependence  of  tnese  two. 
selections.  In  Reference  2-4,  we  outlined  a  scientific  experiment  (not 
in  hardware)  for  the  purpose  of  systematically  investigating  the  impact 
upon  the  performance  of  active  broadband  disturbance  accommodation  of 
various  assumptions  relating  to  the  selection  of  reduced-order  models  and 
of  active  transducers.  This  ex,  eriment  is  not  a  single  design  attempt, 
but  rather  a  process  involving  a  sequence  of  many  variously-constrained 
designs,  in  order  to  converge  on  an  appropriate  overall  algorithm  for 
broadband  disturbance  accommodation  with  large  flexible  structures. 

Sections  2  and  3  of  the  present  report  constitute  the  initial 
report  of  progress  in  the  conduct  of  the  aforementioned  experiment. 
Section  2  focuses  principally  on  the  theory  and  application  of  the  active 
transducer  selection  strategy  being  employed,  while  Section  3  focuses  on 
the  related  matters  of  controller  design  and  evaluation.  Two  very 
encouraging  trends  have  emerged  from  the  results  reported  herein: 

(1)  Under  appropriate  assumptions,  substantial  reductions  in  the 
number  of  active  transducer  pairs  below  the  number  of  controlled  modes 
are  feasible,  which  implies  substantial  design  flexibility  may  be 
available?  and 

(2)  Using  an  expanding  family  of  evaluation  models,  stability  of 
reduced-order  controllers  over  a  frequency  range  well  beyond  the  range  of 
the  design  model  has  been  demonstrated. 

The  significance  of  item  (2)  is  that  sophisticated  adjustments  of  the 
controller  design  have  been  deliberately  postponed?  this  result  reflects 
principally  the  reward  of  careful  attention  to  active  transducer 
selection. 
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The  remainder  of  Section  2  is  briefly  summarized.  The  active 
transducer  selection  strategy  being  employed  in  the  scientific  experiment 
was  briefly  outlined  in  a  previous  report  [Ref.  2-5,  Sec.  3] .  In  Section 
2.2,  a  complete  account  of  the  theory  that  underlies  this  strategy  is 
given.  In  addition,  all  phases  of  the  strategy  are  fully  discussed  in 
the  light  of  the  theoretical  building  blocks.  In  Section  2,3,  a  number 
of  applications  of  the  selection  strategy  to  ACOSS  Model  No.  2, .under 
varying  assumptions,  is  presented.  In  Section  2.4,  some  principal 
observations  on  the  results  are  made.  An  appendix,  Section  2.5,  contains 
the  proof  of  some  technical  matters  that  arise  in  Section  2.2. 

2.2  Theory:  A  Selection  Algorithm 

Most  strategies  for  active  transducer  selection  reported  in  the 
literature  use  appropriately-defined  indices  of  controllability  and 
observability  to  evaluate  alternatives.  The  selection  strategy  proposed 
and  briefly  outlined  previously  [Ref.  2-5,  Sec.  3]  approaches  the 
selection  problem  from  a  different  direction.  Linear  combinations  of  the 
modal  influence  vectors  associated  with  each  actuator  are  sought  which 
best  approximate  the  modal  coefficient  vectors  of  variables  to  be 
regulated  (e.g. ,  line-of-sight  errors  in  an  optical  system).  The 
principal  theoretical  building  blocks  of  this  strategy  are  the  special 
properties  of  full-rank  matrix  factorizations  (in  particular,  the 
QR-decomposition)  and  the  representation  theory  for  solutions  of  least 
squares  approximation  problems.  In  what  follows,  a  complete  though 
concise  account  of  this  underlying  theory  is  necessary  to  enable  a  proper 
description  of  the  construction  and  operation  of  the  selection  algorithm. 


2.2.1  Fundamental  Concepts 


As  will  be  seen  shortly,  least-squares  approximation  problems  lead 
naturally  and  directly  to  the  Moore-Penrose  generalized  inverse.  This 
construction  exists  for  arbitrary  rectangular  matrices  A:  nxm,  but  is  not 
quite  so  direct  in  the  case  that  the  inequality 

k  =  rank(A)  <_  min  {n,m}  (2-1) 

becomes  strict,.  In  that  case,  the  ability  to  decompose  A  into  factors  is 
extremely  helpful.  A  decomposition  of  the  form 

A  =  BC  (2-2) 


where  E:  nxk  and  Cs  kxm  satisfy 

rank(B)  =  k  =  rank(C) 

with  k  defined  by  Eq.  (2-1),  is  called  a  full-rank  factorization  of  A. 
It  is  well  known  that  every  matrix  has  such  a  factorization  [Ref.  2-6]  . 
Although  such  factorizations  are  not  unique,  the  class  of  all  such 
factorizations  is  generated  from  a  known  factorization  by  the  relations 


B1  =»  BY  ,  C1  *  Y_1C 

for  some  nonsingular  matrix  Y  [Ref.  2-7].  Full-rank  factorizations 
assume  great  importance  in  numerical  linear  algebra  [Ref.  2-8], 

Essential  to  the  selection  algorithm  to  be  described  later  is  the  0R- 
decomposition  [Ref.  2-9],  which  always  has  embedded  within  it  a  full-rank 
factorization. 
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2 . 2 . 1 . 1  Pile  QR-Decom position 


The  classical  Gram-Schmidt  orthugonalization  [Ref.  2-10,  pp. 
127-126]  is  a  constructive  process  of  transforming  a  set  of  linearly 
independent  vectors  in  an  inner-product  space  into  an  orthonormal  set. 
The  matrix  representation  of  the  result  is  the  QR-de composition.  It 
follows  that  every  matrix  has  a  QR-decomposition,  whose  specific 
structure  depends  on  the  rank  of  the  matrix.  Since  the  results  are 
thoroughly  discussed  in  the  literature  [e.g..  Ref.  2-11],  the  parts 
relevant  here  are  simply  quoted. 

Theorem  2-1  (Existence  of  QR-Decomposition) .  Assume  that  the  matrix  A: 
nxm  has  rank  k=min{n,m}.  Then  there  exists  a  (nonsingular)  column 
permutation  matrix  P:  mxm,  a  unitary  matrix  Q:  nxn,  and  a  matrix  R:  nxm 
containing  an  upper  triangular  submatrix  R^ :  kxk  with  positive  diagonal 
elements,  such  that 


AP  =>  QR 


The  structure  of  the  R  matrix  is 


CM 


(2-3) 


(2-4a) 


( 2-4b) 


or 


[r,  :  s  ] 


( 2-4c ) 


according  as  A  is  determinate  (k=n=m),  overdetermined  (k=m<n),  or  under¬ 
determined  (k=n<m),  respectively.  □ 
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The  word  "unitary"  refers,  in  the  complex  case,  to  a  matrix  allowing 
complex  entries  whose  columns  are  orthogonal  relative  to  the  complex 
inner  product;  the  corresponding  property  in  the  tesj.  case  refers  to  a 
matrix  restricted  to  real  entries  whose  columns  are  orthogonal  relative 
to  the  real  inner  product  ("orthogonal"). 

,  Theorem  2-1  has  an  obvious  extension  to  the  rank-deficient  case 
that  k  <  min{n,m}.  The  extended  result  is  used  only  in  the  proof  of 
Proposition  2-8  in  the  Appendix  (Section  2.5).  The  decomposition  of  Eq. 
(2-3)  is  a  full-rank  factorization  (cf.  Eq.  (2-2))  except  in  the  over¬ 
determined  case  (cf.  Eq.  (2-4b)).  In  that  case,  Eq.  (2-3)  is  equivalent 
to  the  full-rank  factorization 

AP  »  Q1R1 

where  Q-)  is  the  partition  of  Q  compatible  with  R-|  . 

For  the  intended  application,  stronger  properties  of  the 
nonsingular  submatrix  R-j  are  required  than  are  evident  in  the 
Gram-Schmidt  interpretation  as  reflected  in  Theorem  2-1 .  A  powerful  tool 
for  accomplishing  this  is  the  following  triangularization  lemma  of 
Householder  [Ref.  2-12]. 


Lemma  2-2.  Assume  that  n-vectors  a*0  and  v  are  specified.  Then  there 
exists  a  unit  vector  w:  nxl  such  that 

H(a)  =  (I  -  2  ww*)a  =  lajv  □  (2-5) 

n 

In  Eq.  (2-5),  In  is  the  nxn  identity  matrix,  the  asterisk  (*)  denotes 
the  complex  conjugate  transpose  operation,  and  the  double  bars  denote  the 
vector  norm  generated  by  the  inner  product: 
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(a, a) 


* 

a  a 


las 


2  A 


Successive  application  of  Householder  transformations,  each  having  the 
form  defined  in  Eq.  (2-5),  to  columns  of  the  A  matrix,  including 
appropriate  rearrangement  (pivoting)  of  the  columns  at  each  step,  leads 
to  the  desired  triangularization.  Reliable  and  widely  available  software 
for  computing  a  QR-decomposition  uses  the  Householder  lemma  [Ref.  2—1 3J . 
The  stronger  property  that  results  is  the  following. 


Theorem  2-3.  In  Theorem  2-1,  elements  of  the  matrix  R-j  =  [r^j  ]  may 
be  taken  to  have  the  diagonal  dominance  property: 


>  I 


j 


,  j=£,  ; 


£=1 , •••,k 


(2-6) 


In  particular,  this  implies  the  (weaker)  property  that  the  diagonal 
elements  form  a  nonincreasing  sequence: 


r  „  >  r  „  ,  „  ,  > 

££  “  £+1 ,  Jt+1 


JM  ,  •• »,k-1 


(2-7) 


The  subtle  distinction  between  the  two  theorems  above  should 
be  carefully  noted.  The  results  of  Theorem  2-1  describe  every 
QR-decomposition;  the  sharper  results  of  Theorem  z-3  describe  only  those 
QR-decompositions  generated  by  the  use  of  Householder  transformations. 

2. 2. 1.2  Solutions  to  Least-Squares  Problems 

Least-squares  problems  for  finite-dimensional  linear  systems  of 
the  form 
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Ax 


b 


(2-8) 


with  A:  nxm  become  significant  in  the  overdetermined  case  (n>m)  in  which 
there  are  too  many  equations  for  a  solution  to  exist  (in  general). 
Instead,  one  attempts  to  minimize  a  positive-definite  quadratic  function 
of  the  difference,  such  as 

J(x>  A,  b)  =  (b-Ax) * (b-Ax)  (2-9) 

or,  more  generally  (apparently), 

J  (x;  A,  b)  =  (b-Ax) *W (b-Ax)  (2-10) 

V»  t 

where  W:  nxn  is  Hermitian  (i.e.,  W*=W)  and  positive  definite.  Being 
positive  definite  however,  W  can  be  written  in  the  factored  form 

W  =  L*L 

with  Ls  nxn  nonsingular  (Ref.  2-10,  pp.  139-140].  This  enables  Eq. 

(2-10)  to  be  written  as 

Jw(x;  A,  b)  a  (b  -  Kx)*(b  -  Ax)  =  J(x;  A,  b)  (2-11) 

where  b  «  Lb,  A  =  LA,  so  that  the  least-squares  problems  associated  with 
Eq.  (2-10)  and  Eq.  (2-9),  respectively,  are  in  fact  equivalent  under  a 
simple  scaling.  It  is  therefore  sufficient  to  discuss  the  relevant  parts 
of  the  least-squares  theory  using  Eq*  (2-9),  even  though  the  selection 
algorithm  is  formulated  tsing  Eq,  (2-10).  Only  the  unconstrained 
minimization  of  Eq*  (2-9)  is  discussed.  Full  details,  including 
treatment  of  constrained  minimization,  can  be  found  in  Reference  2-11. 

The  minima  of  Eq,  (2-9)  are  precisely  the  set  of  solutions  to  the 
normal  equations 
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0 


(2-12) 


A*(b-Ax)  = 

It  is  worth  noting  that  this  result  is  a  consequence  of  the  geometric 
structure  of  inner-product  spaces  [Ref.  2-7]  and  does  not  follow  from  the 
methods  of  calculus  (the  operation  x  x  of. complex  conjugation  being 
nowhere  differentiable).  In  addition,  Eq.  (2-12)  provides  a  useful 
geometric  insights  at  any  least-squares  minimum,  the  residual  vector 

A  . 

p  =  b  -  Ax 

associated  with  the  linear  system  Eq.  (2-8)  is  orthogonal  to  the  range 
space  (i.e. ,  the  subspace  spanned  by  the  columns)  of  A,  denoted  by 
col(  A) . 

The  solutions  of  the  normal  equations  form  an  affine  space  of 
dimension  m-k,  where  m  and  k  are  defined  by  Eq.  (2-1).  Written 
explicitly,  they  are: 

x  =  A+b  +  (I  -  A+A)aj  ,  us:  nxl  arbitrary  (2-13) 

m 

where  the  superscript  (t)  represents  the  Moore-Penrose  inverse  [Refs. 
2-14,  2-15].  The  basic  properties  of  this  generalized  inverse  will  be 
used  without  elaboration,  as  they  are  thoroughly  documented  [Ref.  2-16] 
and  have  been  discussed  previously  in  a  related  context  [Refs. 

2-17,  2-18].  It  follows  from  these  properties  that  the  two  terms  in  Eq. 
(2-13)  are  orthogonal  and,  therefore,  that  the  particular  least-squares 
solution 


0 

x 


a  A+b 


has  (uniquely)  the  smallest  norm  of  all  the  solutions  Eq.  (2-13). 
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Each  least-squares  solution  gives  the  minimum  value 


J  .  =  1(1  -  AA+)bJ2  (2-14) 

min  n 

for  the  functional  Eq.  (2-9).  Geometrically,  this  represents  the  length 
(squared)  of  the  projection  of  b  onto  col(A)J\  the  orthogonal  complement 
of  the  range  space  of  A.  In  particular,  it  follows  that  the  minimum 
valuq  is  zero  if  and  only  if 


b  s  col(A)  (2-15) 

i.e.,  that  Eq0  (2-8)  has  an  exact  solution. 

'he  representation  of  least-squares  solutions  can  be  sharpened 
substantially  if  the  Moore-Penrose  inverse  is  displayed  using  the  QR- 
Jecomposition  of  A.  Two  cases  are  distinguished: 

(I)  k  a  n  <_  m,  for  which  Eq.  (2-15)  is  true;  and 
(IT)  k  =  m  <  n,  for  which  Eq.  (2-1 5) ‘is  false. 


Theorem  2-4  (Case  I).  Assume  that  A  has  the  decomposition  of  Eq.  (2-3) 
with  k  =  n  <_  m.  Then  the  minima  of  Eq.  (2-9)  are  given  by: 

(p  ^x)  =  R+Q  b  +  (I  -  R+R)uj  ,  id:  nxl  arbitrary  (2-16) 

If  n=m,  the  set  of  Eq.  (2-16)  collapses  to  the  unique  solution 

(P_1x)  =  R~Vb  •  (2-17) 

In  both  instances,  the  minimum  value  in  Eq.  (2-14)  is  zero.  □ 

This  result  is  used  in  the  initial  stages  of  the  selection  algorithm. 
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an  even  richer  structure  develops  in  the  next  case. 

Theorem  2-5  (Case  II).  Assume  that  A  has  the  decomposition  of  Eq.  (2-3) 
with  k  =  m  <  n.  Denote  by  Q  =  [Q-j  :  Q2]  the  partition  of  Q  compatible 
with  the  partitior  of  R  in  Eq.  (2-4b).  Then  the  minimum  of  Eq.  (2-9)  is 
given  (uniquely)  by: 

(p-’x)  =  r”1  Q1  *b  (2-18) 

and  the  minimum  value  in  Eq.  (2-14)  is 

Jmin  '  ,a2*b'2  °  l2-,9> 

The  significance  of  Theorem  2-5  arises  from  the  riqh  structure  of  the 
mappings  associated  with  the  matrices  Qi :  nxk  and  Q2:  nx(n-k),  which 
are  briefly  summarized  next.  They  all  follow  from  viewing  Eq.  (2-3)  in 
the  equivalent  form 

Q*( AP)  =  R  (2-20) 

which  exhibits  the  final  result  of  the  Householder  construction  (cf. 

Lemma  2-2).  It  follows  immediately  that  Eq.  (2-20)  embodies  the  two 
equations: 

Qn*(AP)  »  R  •  (2-21) 

and 
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0 


(2-22) 


q2  (AP)  = 


From  Eqs.  (2-21)  and  (2-22)  it  follows  quite  directly  that  Q-|  and  Q2 
are  orthogonal: 


*  * 

Q_  Q-  =  0  ,  Q„  Q.  =0  (2-23a) 

*1  *2  2  1 

«t*«,  -  h  '  %‘«2  "  Vk  (2-23b) 


and  have  range  spaces 

col(Q  )  a  col(AP),  col(Q2)  =>  col{  AP)  (2-24) 


respectively.  Furthermore  the  product  mappings 

W  3 

and 

W  3 

are  orthogonal  projections,  respectively  (of.  Eqs.  (2-23),  onto  the  range 
spaces  listed  in  Ea.  (2-24).  Note  that  Eq.  (2-26;  provides  the  link 
between  Eq.  (2-14)  and  the  result  in  Eq.  (2-19).  Theorem  2-5  forms  the 
principal  basis  for  the  selection  algorithm. 

The  behavior  of  the  minimum  value  in  Eq.  (2-14)  as  a  function  of 
A,  and  in  particular  as  columns  are  removed  from  A,  is  of  considerable 
interest  for  the  active  transducer  selection  problem.  Roughly,  one 
expects  that  the  minimum  value  does  not  change  it  columns  of  A  are  simply 
rearranged,  and  that  it  does  not  decrease  if  a  column  of  A  is  removed. 


Xn  -  Q2Q2 


Jn  -  Vi 


(2-25) 


(2-26) 
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Theorem  2-6 .  Let  p  be  an  integer,  1  £  p  <_  m,  and  denote  by 
an  arbitrary  p-column  submatrix  of  A.  Then 

*  P 


111,1  J(x;  A,  h)  <  mln  J(y;  t 
xeEm  yeE^ 


f  l 

\  P 


,  b) 


(2-27) 


wv3re  E°  denotes  the  appropriate  a-dimensional  i  n.tAt  tvwdact  space.  Cl 


The  proof  of  this  result  is  deferred  to  the  Appendix  (Section  2.5). 

2.2.2  Algorithm  Formulation 

It  is  now  possible  to  give  a  precise  formulation  of  the  acrive 
transducer  selection  algorithm.  It  is  assumed  in  what  follows  that: 

(1)  A  reduced-order  structural  model  has  been  chosen;  and 

(2)  Sensors  are  to  be  colocated  with  actuators. 

The  focus  is  on  determining  an  appropriate,  and  relatively  small, 
selection  of  actuators  from  a  prespecified  collection  of  candidates. 

Dynamic  equations  in  modal  coordinates  for  the  reduced-order  model 
to  be  considered  have  the  general  form: 

n  +  2zn.n  s-  f)2n  =  (*Iba)ua  +  ($tbd)ud  (2-28) 

/  T  T>t  /  T  T\T  • 

y  -  (*  v)  n  +  (♦  V)  n  (2-29) 

r  T  T\T 

z  =  (*  D  )  n  (2-30) 
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where  n  =  (hif,*,/hn)T  is  the  vector  of  modal  coordinates  retained 
in  the  reduced-order  model?  ft  =  diagtajp):  nxn  is  the  matrix  of 
characteristic  frequencies;  $:  vxn  is  the  truncation  of  the 
principal-axis  matrix  of  transformation  defining  the  modal  coordinates 
(mode  shapes  of  the  reduced-order  model) ,  v  being  the  number  of  physical 
generalized  coordinates?  Z  =  diag(  ;■),*•  v?n)  :  nxn  is  the  matrix  of 
assorted  modal  damping  ratios  representing  inherent  structural  damping; 
uA  =  (ui/»*»/um)T  is  the  vector  of  inputs  to  the  actuators? 
up  S  (w-i ,  is  the  vector  of  disturbance  inputs? 

y  =  t  is  the  vector  of  outputs  from  the  sensors? 

z  =  (zi,»«*,zg)T  is  the  vector  of  variables  to  be  regulated;  and 
Ba:  vxm,  Bq?  vxy,  Cp:  Jlxv,  Cy?  £xv,  Ds  3x\?  are  the  influence 
matrices  associated  with  the  actuators,  disturbances,  displacement 
sensors,  rate  sensors,  and  regulated  variables,  respectively.  The 
superscript  (T)  denotes  matrix  transpose.  Equations  (2-28)  through 
(2-30)  are  written  in  a  form  that  highlights  the  presence  of  generalized 
node  shapes  (modal  influence  vectors)  of  the  form  $Tx. 

The  selection  problem  is  defined  in  terms  of  two  sets  of 
parameters? 

(1) -  A  collection  of  actuator  candidates,  consisting  of  the 

T  i  i 

generalized  node  shapes  $  b:,  where  b^  is  the  j-th  column  of 

Pi  Pi 

the  matrix  BA  in  Eq.  (2-28)?  and 

(2)  A  collection  of  regulation  target  vectors,  consisting  of 

T  T 

tne  generalized  node  shapes  $  d^,  where  d^  is  the  i-th  row 

of  the  matrix  D  in  Eq.  (2-30). 

It  is  assumed  that  the  designer  has  a  relatively  wide  range  of  choices 
for  actuators  (in  particular,  having  redundancy)  both  by  functional  type 
and  location  for  inclusion  in  the  candidate  collection  .  The  idea 
behind  collection  dR  is  as  follows:  actuators  are  to  be  selected  from 
the  collection  which  have  a  relatively  strong  influence  on  each  of  the 
vectors  in  collection  £8  ,  thereby  enhancing  the  prospects  of  achieving 


2-14 


the  desired  regulation  of  the  vector  z  in  Eq.  (2-30) .  In  order  to 
clarify  ideas  and  simplify  the  notation,  the  remaining  discussion 
identifies  the  regulated  variable  z  with  the  two-vector 
( Elosx *  eL0SY )  that  is  of  interest  in  the  applications'  of  Section  2.3 
and  which  represents  line-of-sight  (LOS)  error  parameters  for  an  optical 
system.  The  collection  .35  then  consists  of  the  two  elements: 


,  LOSX  A 
b  = 


T 

i  d, 


LOSX 


LOSY  A 
b  = 


«  d. 


LOSY 


(2-31) 


T  T 

where  d  and  d  are  the  rows  of  D  corresponding  to  the  regulated 
LOSX  LOSY 

variables  Elosx  and  eL0SY  EcI*  (2-30). 

Using  the  parameters  from  collections  *.*/  and  &  ,  an  ideal 
selection  problem,  and  a  least-squares  approximation  to  it,  are  defined. 
Denote  by 


A  =  [jrf ]  =  [a1  •  •••  |  a*1]  :  n  x  p  (2-32) 

the  matrix  formed  from  elements  of  the  actuator  candidate  collection  jl . 
(The  type  of  correspondence  defined  by  Eq.  (2-32)  is  used  repeatedly  in 
what  follows.)  It  is  assumed  that  the  collection  jd  is  rich  enough  to 
span  the  space  En  of  modal  coordinates;  i.e. , 


sp(^)  =  col(A)  =  En 


(2-33) 


Then  there  exist  coefficient  vectors  x  =  (x 
such  that 


• • • ,x  ) 

ft 


y  =  (y 


r 


^LOSX 


Ax 


«  l 


.  ,  a^x. 
3=1  3 


and 
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LOSY 

b 


Ay 


Ideal  Selection  Problem:  Find  a  minimal  subcollection  c*/'  of  and 
real  coefficient  vectors  x*  =  (Xj),  y*  5  ( y j ) »  j  e  lijd'),  such  that 


LOSX 


=  A'x'  =  y.  .  a-’x'  . 

TSl(  jl'  )  3 


and 


LOSY  . .  ,  _  v  j  , 

5  “  A  y  ■  ^jelt  jl' )  a  r) 


where  the  index  set  H.jd')  is  the  subset  of  { 1 , * • • # y }  that  identifies 
which  elements  of  <jd  are  retained  in  -J&' ,  and  A'  = 

as  a  more  feasible  alternative  to  the  ideal  selection  problem,  the 
following  least-squares  approximation  is  formulated.. 


Approximate  Selection  Problem;  Find  a  subcollection  of  and 
corresponding  coefficient  vectors  x'  =  (Xj),  y'  3  ( y  j ) ,  j  e  I(^') 
such  that 


...  ,  .  .  -  .LOSX, 

(1)  Jw(x';  a',  b  ) 


,,  LOSX  ,  .  .  <  T  ..LOSX 

(b  -  A'x' )  W (b 


A'x')  and  (2-34) 


_  .  ,  At  .LOSY,  A 

Jw(y';  a',  b  )  = 


..LOSY  _  .  .  ,T„,,LOSY 
(b  -  A'y' )  W(b 


A'y'  ) 


(2-35) 


are  minimized,  where  A'  =  and 

(2)  the  subcollection  v*/1  is  minimal  with  respect  to  acceptable 
increases  in  the  minimum  values: 


J  .  (A', 
min 

LOSX . 

D  ) 

A 

min 

X’ 

V*'! 

A', 

bL°SX) 

(2-36) 

J  .  (A', 
min 

bL°SY) 

A 

min 

y' 

yy. 

A’, 

bL0SY) 

(2-37) 

In  the  above,  the  matrix  W:  nxn  is  chosen  to  be  real,  symmetric,  and 
positive  definite.  However,  in  view  of  the  equivalence,  under 
appropriate  scaling,  conveyed  by  Eq.  (2-11),  the  discussion  below  omits 
explicit  reference  to  the  weighting  matrix  W.  Note  that,  although  the 
correspondence 


ot'  +  [jtf'l  =  A' 

defined  by  Eq.  (2-32)  does  not  distinguish  between  distinct  column 
rearrangements,  Proposition  2-7  (Section  2.5)  shows  that  such 
distinguishability  is  inessential  for  its  intended  usage.  Thus,  the 
approximate  selection  problem  is  well-defined. 

A  brief  overview  of  the  algorithm  for  the  solution  of  the 
approximate  selection  problem  is  as  follows.  A  sequence  of  eliminations 
from  the  actuator  candidate  collection  is  generated  on  the  basis  of 
information  from  least-squares  solutions  corresponding  to  Bqs.  (2-34)  and 
(2-35)  at  each  step.  These  solutions  are  displayed  using  the  appropriate 
QR-decomposition.  Note  that  each  elimination  also  leads  to  a  reduction 
in  dimension  of  the  inner-product  space  in  which  the  least-squares 
minimization  is  done.  No  matter  what  elimination  sequence  develops, 
Theorem  2-6  guarantees  a  sequence  of  least-squares  approximations  that  is 
steadily  deteriorating  in  quality  in  the  sense  that  the  sequence  of  norms 
for  the  residual  vector  at  each  least-squares  solution  is  nondecreasing. 
When  has  oeer.  reduced  by  elimination  to  a  collection  for  which 

further  elimination  (say  to  o^")  leads  to  the  condition  that  either  of 
the  minimum  values  Eq.  (2-36)  or  (2-37)  (corresponding  to  <.*/"  rather  than 
to  )  exceeds  a  specified  tolerance  tq»  the  elimination  sequence  is 
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terminated,  and  is  taken  to  be  a  minimal  subcollection.  Next, 
details  of  the  elimination  sequence  are  outlined. 

2. 2. 2.1  Criteria  for  Elimination 

In  view  of  Theorem  2-4,  the  search  for  a  minimal  subco} lection  in 
the  sense  of  the  approximate  selection  problem  eventually  requires 
criteria  for  intelligently  eliminating  a  column  from  a  determinate  or 
overdetermined  matrix.  Two  such  criteria  are  discussed  prior  to 
describing  the  elimination  sequence.  Consider  a  stage  of  the  sequence  at 
which  the  matrix  A  =  of  the  remaining  actuator  candidates  is 

(say)  overdeternuned  with  £  <  n  columns.  The  QR-decomposition  of  A  has 
the  form  of  Eq.  (2-3)  with  R  having  the  form  of  Eq.  (2-4b).  Denote  by  x 
and  y,  respectively  the  unique  solutions  (cf.  Eq.  (2-18))  to  the 
least-squares  problems  associated  with  Bqs.  (2-34)  and  (2-35).  The 
criteria  for  elimination  of  actuator  candidates  from  A  are  defined  in 
terms  of  the  matrices  R,  x,  and  yj  in  particular: 

(I)  the  relative  magnitude  of  the  smallest  diagonal  element 
of  Ri ;  and 

(II)  the  relative  magnitude  of  the  products  r^x^,  i*1, •••,£, 
andrjjyj,  j=1,  •••,£,  respectively. 

The  rationale  for  each  criterion  is  briefly  discussed. 

Criterion  I.  The  scalar  equations  for  the  QR-decomposition  of  A  have  the 
triangular  structure: 


ai  -  “,r,j  +  <!%j +  —  **Sj  ■  >’••••'* 


(2-38) 


where  the  diagonal  elements  r j j  of  R^  form  a  nonincreasing  sequence 

(cf.  Eq.  (2-7)).  Denote  by  S£_i  the  ( £-1  )-dimensional  subspace  spanned 

1  £-1  £  £ 
by  the  vectors  {q  ,•*•,  q  }.  Since  the  vector  a  -q  r  is  m  5  , 

it  *4  £**  * 

r  is  an  upper  bound  for  the  Euclidean  distance  between  a.  and  S „ 

£2  £  £-1 
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d(a  ,  S^_1  )  =  inf  { Da  — s  II :  s  £  S^_1  } 


(2-39 


< 


'U 


Because  of  the  triangular  structure  of  Eas.  (2-38),  S^_i  is  also  the 
linear  span  of  {a1,***,  a*-1 }.  Therefore,  Eq.  (2-39)  is  a  measure  of 
the  near-dependence  of  a&  upon  the  preceding  column  vectors  a^, 
i=1  , « . £-1 .  Sufficiently  small  values  of  r%£  relative  to  r£_if£_i 
suggest  that  might  be  eliminated  without  an  unacceptable  increase  in 
the  resulting  minimum  value  Eq.  (2-36)  and  (2-37)  associated  with  the 
truncated  matrix. 


Criterion  II.  Eliminations  in  accordance  with  Criterion  I  are  made  only 
on  the  basis  of  the  QR-decomposition  of  A,  and  without  regard  for  the 
least-squares  solutions  that  determine  the  quality  of  the  desired 
approximation.  Criterion  II  takes  those  solutions  into  account.  In 
comparing  the  numbers  |  r  |  »  i=\, ••♦,£,  suppose  that  the  minimum 
occurs  at  some  index  o  <  .It  i.e.. 


r  x 
00  0 


min 

i=1 


(2-40: 


Recalling  the  nonincreasing  property  Eq.  (2-7)  of  the  diagonal  elements 
rj_jL,  it  follows  from  Eq.  (2-40)  that 


min 

i=o, •••,£ 


( 2  -  4 1 : 


Equation  (2-41)  is  an  indication  that  column  a 0  is  less  important  to 
the  least-squares  approximation  than  the  succeeding  columns  a3, 
j=o+i, Sufficiently  small  values  of  the  minimum  in  Eq.  (2-40) 
relative  to  the  other  |rj^Xj.|  values  suggest  that  a0  might  be 
eliminated  from  A.  Similar  considerations  relating  to  the  numbers 
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4 


f 


fi 


|rj j^j I  '  3=1 '***'£'  suggest  the  elimination  of  a  (possibly  different) 
column  aY  from  A.  Taken  together,  Criteria  I  and  II  may  suggest  up  to 
three  distinct  candidates  for  elimination  from  A. 

2. 2. 2. 2  The  Elimination  Sequence 

Having  specified  the  following  parameters: 

(1)  an  actuator  candidate  collection  jd ; 

(2)  a  regulation  target  collection  3B  (e.g. ,  {bL0SX,  bL0SY}); 

(3)  a  least-squares  weighting  matrix  W  (notational  references 
suppressed);  and 

(4)  a  least-squares  minimum  tolerance  parameter  xq  >  0; 
the  elimination  sequence  can  be  defined. 


Step  1 :  Initial  Factorization.  Set  A^  =  A  =  \jl  1  and  compute  a  QR- 
decomposition  of  A^0^ .  [Because  of  the  assumption  Eq.  (2-33),  the 
R-matrix  in  the  decomposition  (cf.  Eq.  (2-3)) 


,  Qoyi> 


has  the  structure  of  Eq.  (2-4c).J 


(2-42) 


Step  2:  Initial  Least-Squares  Solutions.  Denote  A^  =  A^P^  and 
solve  the  least-squares  problems  of  the  form  Eqs.  (2-34)  and  (2-35), 
respectively,  for  A^1).  [It  follows  from  Theorem  2-4  that  there  are 
infinitely  many  least-squares  solutions  (cf.  Eq.  (2-16)),  but  all  give 
zero  for  the  minimum  values  Eqs.  (2-36)  and  (2-37)  corresponding  to 
A<1>.] 

Step  3:  Initial  Elimination.  Truncate  A^1)  by  retaining  only  its 
first  n  columns;  denote  the  result  by  A  ( 2 )  •  i.e., 


(2) 


=  ( 


,Pl<n) 


n  x  n 
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where  pi  is  the  permutation  mapping  of  the  indices  {l,»«»,u}  induced  by 
the  column  permutation  matrix  p(')  of  Step  2.  [Because  of  the 
assunption  Eq.  (2-33)  and  the  structure  of  the  matrix  R^  ) ,  the  first  n 
columns  of  A^ 1  ^  are  linearly  independent,  Moreover,  A^)  inherits 
the  QR-decomposition 


,  (2) 


=  a'"*!” 


from  Eq.  (2-42)  in  which  the  R-matrix  has  the  form  of  Eq.  (2-4a)'.  There 
are  unique  solutions  x(2)  and  y(2)  to  the  least-squares  problems  of 
the  form  Eqs.  (2-34)  and  (2-35),  respectively,  for  A^)  (cf.  Eq. 

(2-17)),  and  the  corresponding  minimum  values  Eqs.  (2-36)  and  (2-37)  are 
again  zero.] 


Step  4:  Secondary  Elimination.  Examine  the  numbers  r^j  r!Vxf2\ 

- e. -  ,  nn  li  i 

i*1,»»»,n»  and  r^y^,  j=l,.«.,n,  associated  with  A^,  x^2\  and  y*2\ 

Using  either  Criterion  I  or  II,  select  a  column,  say  the  c-th,  of  A^) 

for  elimination.  Truncate  A^2^  by  removing  the  a-th  column.  [When  the 

criteria  suggest  more  than  one  column  of  A^2^  for  elimination  ,  all 

eliminations  may  be  done  in  parallel.  The  result  is  a  (generally 

expanding)  lattice  of  alternatives  for  elimination.] 


Step  5;  Secondary  Factorization.  Denote 


(3) 


4  [.?!<’ 


•  •  •  •  a 


Pl(0)-1  ;  Pl(o)+1 


]s  (n-1) 


x  n 


and  compute  a  QR-decomposition  of  A^2^.  [The  columns  of  A^2^  are 
linearly  independent,  having  been  inherited  from  a(2),  The  QR- 
decomposition 


2-2] 


S(3V2> 


Ql2V2) 


has  an  R-matrix  with  the  structure  of  Eq.  U-'lb),] 

Step  6:  Secondary  Least- Squares  Solutions.  Denote  A^  »  &^)p(2)  and 
solve  the  least-squares  problems  of  the  form  Eqs.  (2-34)  and  (2-35), 
respectively,  for  A<4).  [It  follows  from  Theorem  2-5  that  each  of  the 
two  least-squares  problems  has  a  unique  solution  (cf.  Eq.  (2-18)),  and 
the  minimum  values  corresponding  to  Eqs.  (2-36)  and  (2-37),  respectively, 
are  generally_jionzero  (cf.  Eq.  (2-19)).) 

Step  7:  Subsequent  Eliminations.  Repeat  Steps  4  through  6, 
successively,  applying  Step  4  of  the  current  cycle  to  the  results  of  Step 
6  for  the  preceding  cycle. 

Step  8:  Termination.  Assume  that  the  preceding  eliminations  have 
reduced  the  initial  collection  to  the  subcollection  and  that  each 

elimination  (say  to  jd" )  from  jl'  via  the  application  of  Steps  4  through 
6  leads  to  corresponding  minimum  values  (cf.  Eqs.  (2-36),  (2-37))  such 
that  either 


or 


J_in  (A"'  bL°SX)  >  To  (2-43) 

IiOSY,  ... 

J  (A",  b  )  >  T_  (2-44) 

min  0 


where  A"  =  [**/"].  Terminate  the  elimination  sequence;  ud'  is  a  minimal 
selection. 

The  actual  numerical  implementation  of  the  termination  tests  Eqs. 
(2-43)  and  (2-44)  is  modified  slightly  when  very  small  residual  vectors 
in  the  least-squares  solutions  are  required,  as  in  the  applications  -o  be 
presented  subsequently.  The  test  Eq.  (2-43)  has  the  form 
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(2-45) 


r  n  2 

U=1  pi  >  T0. 

where  p  =  (pi,»**,pn)T  is  the  residual  vector  bL0SX-A''x"  at  a 
least-squares  solution  x"  corresponding  to  A".  The  actual  computation  of 
the  squares  of  small  numbers  in  Eq.  (2-45)  in  undesirable.  In  view  of 
the  equivalence  of  all  norms  on  En,  and  in  particular,  from  Jensen's  - 
inequality  {Ref.  2-19]: 


max  .|  i  '  /rn  i  i2J/2  -  n  ,  , 

i»l ,  • »»,n  |pit  •<-  ^i-1  |pi|  )  1  Ei-i  |pi|  1  nM 


the  test  Eq.  (2-45)  is  equivalent  to  the  condition 


i=>1 , •  •  • ,n 


(2-46) 


for  appropriately  chosen  T|.  Similar  remarks  apply  to  the  test  Eq. 
(2-44).  In  addition  to  being  simpler  and  subject  to  less  computational 
error,  Eq.  (2-46)  also  provides  information  (not  provided  by  Eq.  (2-45)) 
as  to  which  of  the  modal  coordinates  is  associated  with  the  unacceptable 
approximation. 

It  should  be  observed  that  the  elimination  sequence  has  embedded 
within  it  a  considerable  degree  of  design  flexibility  that  goes  far 
beyond  the  choice  of  the  parameters  jl,  SB  ,  W ,  and  to*  Step  3 
accomplishes  a  (generally  quite  large)  reduction  in  the  number  of 
actuator  candidates  from  y  (cf.  Eq.  (2-32))  to  n.  The  candidates  that 
are  retained  in  Step  3  are  determined  by  the  column  permutation  pt1) 
generated  in  Step  1.  That  permutation  in  turn  is  the  result  of 
successive  applications  of  the  Householder  Lemma  2-2  that  are  performed 
in  Step  1  in  order  to  achieve  the  diagonal  dominance  property  Eq.  (2-6) 
in  the  factorization  of  Eq.  (2-42).  In  general,  p(')  is  not  unique  in 
this  regard.  The  implications  of  this  lack  of  uniqueness  have  yet  to  be 
fully  explored. 
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Applications  with  ACOSS  Model  No.  2 


In  the  description  of  the  scientific  experiment  for  active  control 
synthesis  with  generic  structural  models  [Ref.  2-4),  four  principal 
elements  of  the  synthesis  process  are  identified  as  variables  of  the 
experiment: 

(1)  Design  of  the  basic  structure  to  be  controlled; 

(2)  Selection  of  reduced-order  models; 

(3)  Selection  of  active  transducers;  and 

(4)  Determination  of  the  controller  feedback  structure. 

For  the  work  summarized  in  the  present  report,  the  first  and  last  of 
these  elements  are  fixed  in  order  to  focus  on  the  other  two,  especially 
active  transducer  selection. 

The  first  element,  basic  structural  design,  is  a  significant 
variable  in  relation  to  the  generic  structural  models  because  of  the  many 
modifications  to  original  designs.  In  particular,  the  original 
structural  design  for  ACOSS  Model  No.  2  [Ref.  2-20)  has  undergone  three 
major  revisions  [Ref.  2-21).  The  choice  of  Revision  1  for'  the  work 
reported  here  represents  a  compromise  which  retains  improvements  on  the 
original  design  in  the  modeling  of  the  attachments  for  certain  rigid- 
body  elements,  while  avoiding  the  extremes  of  very  high  inherent 
flexibility  (Revision  3)  or  very  high  inherent  stiffness  (Revision  4). 

The  last  element,  controller  feedback  structure  determination,  is 
held  fixed  for  the  work  reported  here  in  order  to  avoid  obscuring  th» 
view  of  whatever  improvements  in  performance  may  be  attributable  to 
judicious  selections  of  reduced-order  models  and  of  active  transducers. 
The  choice  of  the  disturbance-rejection  control  feedback  structure  [Re.', 
2-22;  Sec.  3.6),  which  assumes  full  knowledge  of  the  statistics  of  the 
disturbance,  provides  an  upper  bound  for  control  system  performance. 

The  disturbance  to  be  accommodated  is  assumed  to  be  an  aperiodic 
disturbance  with  a  constant  power  spectral  density  at  frequencies  from  0 
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to  5  Hz,  applied  in  a  fixed  direction  at  nodes  no.  37  ard  46  of  the 
structural  model.  Full  details  of  the  disturbance  definition  and 
relevant  assumptions  have  been  given  previously  [Ref.  2-5,  Sec.  2). 

To  assist  in  the  selection  of  reduced-order  models,  criteria  for 
ranking  the  structural  modes  are  commonly  used  [e.g. ,  Ref.  2-2].  In  the 
present  context,  ranking  of  modes  is  done  by  comparing  open-loop  modal 
responses  of  the  regulated  variable  to  the  disturbance  input;  i.e.,  the 
root-mean-square  4-RMS)  value  of  the  response  of  the  two-input,  single¬ 
output  system  between  the  disturbance  input  and.  the  contribution  of  each 
mode  to  the  line-of-sight  error  (cf.  Eqs.  (2-28)  and  (2-30)).  The 
precise  formulation  has  been  given  previously  [Ref.  2-5,  Sec.  5] .  The 
resultant  ranking  of  flexible-body  modes  for  the  chosen  structural  model 
(Revision  1)  is  shown  in  Table  2-1.  An  inevitable  feature  of  any  such 
ranking  procedure  is  that,  in  general,  a  selected  group  G  of 
consecutively-ranked  modes  is  not  contiguous  in  frequency;  i.e.,  there 
exist  modes  not  in  G  which  are  interlaced  in  frequency  with  modes  in  G. 
For  example,  the  first  eleven  modes  according  to  the  ranking  of  Table  2-1 
are  interlaced  in  frequency  with  seven  modes  (11,  15,  17-21)  of 
considerably  lower  rank.  The  interlacing  can  be  seen  clearly  in  the 
frequency  listing  of  Table  2-2.  A  reduced-order  model  having  this 
property  is  said  to  be  .interlaced;  otherwise  i't  is  called  contiguous. 

As  indicated  in  the  preceding  section,  the  objective  of  the  active 
transducer  selection  algorithm  (roughly  speaking)  is  to  choose  a 
relatively  small  number  of  those  actuators  that  are  expected  to  have 
relatively  high  effectiveness  in  achieving  the  desired  regulation  of 
line-of-sight  error.  It  therefore  seems  desirable  that  the  selection 
process  lead  to  a  rich  lattice  of  alternatives  (Step  4)  for  elimination 
from  the  initial  class  of  actuator  candidates.  This  qualitative  property 
is  termed  richness .  It  is  certainly  desirable  that  upon  termination  of 
the  elimination  sequence,  the  number  of  retained  actuators  is 
substantially  smaller  than  the  number  of  modes  to  be  controlled. 

Denoting  by  p  the  number  of  actuators  retained  in  a  minimal  selection, 
the  level  of  reduction 
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(2-47) 


A  =  A  (  js/  »  tq)  =  n-p 

is  an  important  quantitative  measure  of  the  effectiveness  of  the 
selection  algorithm. '  It  is  obvious  that  the  level  of  reduction  should 
depend  on  the  tolerance  parameter  xq',  it  also  turns  out  to  be. a  very 
strong  function  of  the  class  jd  of  actuator  candidates  chosen. 

Principal  observations  with  respect  to  the  performance  of  the 
selection  algorithm  that  have  emerged  from  work  to  date  are  the 
following: 

(1)  Level  of  reduction  increases  with  augmented  richness; 

(2)  Richness  is  augmented  when: 

(a)  Interlaced  reduced-order-design  models  for  which  active 
transducer  selection  is  tailored  are  made  contiguous; 

(b)  The  ciass  of  actua:or  candidates  is  enlarged; 

(c)  Suitable  restrictions  are  imposed  upon  a  sufficiently 
large  actuator  candidate  class; 

(3)  Minimal  selections  which  realize  the  same  level  of  reduction 
are  not  necessarily  unique: 

Four  design  examples  with  ACOSS  Model  No.  2  are  described  next. 
These  examples  demonstrate  the  observations  noted  and  show  in  addition 
some  of  the  design  flexibility  that  is  available  when  relatively  high 
levels  of  redaction  are  achievable. 

2.3.1  Specifications  for  Design  Examples 

Considerations  which  apply  to  each  of  the  four  examples  are  the 
following: 

(1)  The  basic  structural  design  is  Revision  1  of  ACOSS  Model 
No.  2  (Ref.  2-21]  with  Z  =  diag(O.OOl)  in  Eq.  (2-28); 

(2)  A  broadband  disturbance  [Ref.  2-5,  Sec.  2)  is  applied  to 
nodes  no,  37  and  46  of  the  structvure,  assuming  constant 
direction; 
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(3)  Selection  of'  modes  for  inclusion  in  reduced-order  models  for 
controller  design  is  assisted  by  the  disturbance  response 
ranking  of  Table  2-1 ; 

(4)  The  regulation  target  collection  0}  consists  of  the  elements 
defined- by  Eq.  (2-31)  (whose  dimension  is  that  of  the  design 
mode 1 ) ; 

(5)  Ihe  least-squares  weighting  matrix  W  has  the  form: 


W  =  diag(w^ ) :  n  x  n 


(2-48) 


(6) 


where  wj_  is  the  RMS  LOS-error  disturbance  response 
associated  with  mode  m  of  the  reduced-order  design 
model  as  listed  in  Table  2-1 ; 

Ihe  least-squares  minimum  tolerance  parameter  is 


10 


-6 


(2-49) 


The  four  examples  are  distinguished  by  differing  specifications  on 
the  reduced-order  design  model  and  actuator  candidate  class  .  The 
specifications  are  summarized  in  Table  2-3  and  are  discussed  briefly 
below. 


Reduced-order  design  model.  In  Example  A,  the  reduced-order  design  model 
consists  only  of  the  modes  ranked  1  through  11  in  Table  2-1;  it  is 
interlaced.  In  all  the  succeeding  examples,  this  design  model  is  made 
contiguous  Ly  incorporating  the  formerly  excluded  modes  between  mode  7 
and  mode  24.  However,  mode  19  —  a  solar  panel  torsion  mode  —  is 
excluded  from  the  numerical  implementation.  Ihe  associated  entry  in  the 
least-squares  weighting  matrix  W  (cf.  Eq.  (2-48) )  is  extremely  small 
(10"22) 

causing  the  resulting  weighting  matrix  to  be  poorly  conditioned 
numerically  with  respect  to  inversion.  Ihe  small  weighting  is  a 
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reflection  of  the  very  low  influence  of  this  mode  (ranking  at  most  147 
out  ofB  156)  on  any  of  the  modeled  node-connecting  elements. 

Actuator  candidate  class  ,?l  .  In  Examples  A  and  B,  only  axial  member 
actuators  (137  total)  on  node-connecting  elements  of  the  finite-element 
model  are  allowed.  In  Example  C,  the  class  of  candidates  is  expanded  to 
include  translation  actuators  at  all  nodes,  an-aggregate  of  311.  In 
Example  D,  certain  nodal  translation  actuators  are  disallowed;  in 
particular,  these  are  the  z-translation  actuators  on  equipment  section 
nodes  (6  total)  and  both  y-translation  and  z-translation  actuators  on 
interior  solar  panel  nodes  (16  total). 

2.3.2  Selection  Results 

Active  transducer  selection  results  with  Examples  A  through  D  are 
sunmarized  in  Table  2-J;  visual  representations  are  given  in  Figures  2-7 
through  2-9.  The  results  are  discussed  briefly. 

Example  A  yields  somewhat  unsatisfactory  results  from  the 
viewpoint  of  design  flexibility  in  actuator  selection  (A=1 ) •  Only  two 
alternatives  are  available  (Step  4)  for  removing  an  actuator  from  the 
candidates  remaining  after  completion  of  the  initial  elimination  (Step 
3) .  Both  alternatives  lead  to  least-squares  solutions  which  exceed  the 
specified  termination  tolerance  (Eq.  (2-49)).  Hie  least  undesirable  of 
the  two  elimination  alternatives  leads  to  the  minimal  selection  shown  in 
Fig.  2-1.  It  can  be  seen  that  the  measure  of  the  residual  vector  in  Eq. 
(2-46)  slightly  exceeds  the  specified  tolerance  to  in  the  regulated 
varible  associated  with  bL0SX  (cf.  Eq.  (2-31)).  inclusion  of  the 
member  connecting  one  solar  panel  (at  node  52)  to  the  equipment  section 
(at  node  43)  is  somewhat  surprising,  since  in  a  static  comparison  of 
influence  upon  each  of  the  modes  in  the  design  model  (intra-rcw 
comparison  of  entries  in  the  matrix  of  Eq.  (2-28)),  this  member 

ranks  no  higher  than  91st  of  all  137  members. 
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Example  B  yields  a  somewhat  richer  lattice  of  elimination 
alternatives  and  a  substantial  improvement  of  the  level  of  reduction  in 
the  minimal  selection.  This  is  attributable  entirely  to  having 
(minimally)  augmented  the  reduced-order  design  model  to  make  it 
contiguous,  as  Table  2-3  indicates.  The  minimal  selection  is  shewn  in 
Fig.  2—2. 

Example  C  yield,  only  a  slightly  richer  lattice  of  elimination 
alternatives,  but  a  quite  substantial  improvement  in  the  level  of 
reduction  in  the  minimal  selection,  in  comparison  to  example  B.  This  is 
attributable  to  having  enlarged  the  class  of  actuator  candidates  to 
include  translation  actuators  at  nodes,  as  Table  2-3  indicates.  This 
explains  the  presence  of  in  Eq.  (2-47).  The  minimal  selection  is 
shown  in  Figure  2-3.  It  is  noteworthy,  but  not  surprising,  that  no  axial 
actuators  appear  in  the  minimal  selection  even  though  the  candidate  class 
includes  them  all.  In  fact,  all  axial  candidates  are  removed  in  the 
initial  elimination  (Step  3). 

Example  D  provides  somewhat  surprising  and  quite  satisfying 
results.  The  lattice  of  elimination  alternatives  is  extremely  rich  and 
leads  to  an  even  more  substantial  improvement  in  the  level  of  reduction 
in  minimal  selections,  of  which  there  are  two  in  this  case.  A  portion  of 
the  branch  of  the  elimination  lattice  that  leads  to  the  minimal 
selections  is  shown  in  Fig.  2-4.  The  label  "TX-34"  indicates  the  removal 
of  the  translational  (T)  actuation  in  the  x-coordinate  (X)  body-axis 
direction  at  node  34,  other  labels  having  an  analogous  interpretation. 
Parentheses  indicate  eliminations  along  a  portion  of  the  branch  that  do 
not  lead  to  one  of  the  minimal  selections.  Preselections  D-17,  D-14,  and 
D-1 1  are  shown  in  Figs.  2-5  through  2-7,  respectively,  and  minimal 
selections  D-9A  and  D-9B  are  shown  in  Figs.  2-8  and  2-9,  respectively. 
Taken  together,  Figs.  2-4  through  2-9  show  quite  clearly  how  the 
elimination  sequence  (beginning  with  Step  4)  leads  systematically  to 
minimal  selections.  The  steady  increase  in  the  measure  of  Eq.  (2-46)  for 
the  residual  vectors  associated  with  bL0SX  and  bL0SY  (Eq.  (2-31)), 
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respectively,  as  the  elimination  sequence  proceeds  may  be  noted, 
consistent  with  the  expectation  of  Theorem  2-6  (Eq.  (2-27)).  In  this 
example,  the  number  of  actuators  in  a  minimal  selection  is  only  half  the 
number  of  controlled  modes.  This  implies  a  very  desirable  degree  of 
design  flexibility  is  available  (and  may  be  beneficial)  by  employing 
intermediate  (ancestral)  selections  such  as  D-14  or  D-1 1  that  lie  along 
the  p'ath  to  a  minimal  selection.'  h  full  explanation  of  how  the 
restriction  imposed  in  proceeding  to  example  D  leads  to  the  improvements 
observed  is  net  yet  available.  However  it  does  provide  a  substantiation 
of  the  design  flexibility  embedded  within  the  elimination  sequence  that 
is  mentioned  at  the  end  of  Section  2.2. 

Minimal  selections  from  Examples  B,  C,  and  L  are  employed  in  the 
design  oJ:  optimal  linear-quadratic  disturbance-rejection  controllers  as 
reported  subsequently  in  Section  3.  Each  of  the  designs  is  evaluated  in 
the  presence  of  a  family  of  contiguous  reduced-order  evaluation  models 
(cf.  Table  2-2)  with  increasing  ranges  of  frequency.-  Stability  well 
beyond  the  frequency  range  (5  Hz)  of  the  design  model  is  demonstrated  in 
each  case. 

2.4  Summary 

Work  accomplished  to  date  in  connection  with  the  scientific 
experiment  for  active  control  synthesis  with  broadband  disturbances  has 
led  to  substantial  insight  into  a  relatively  new  and  systematic  approach 
to  active  transducer  selection.  This  process  leads  to  selections  which 
at  the  same  time  have  relatively  high  influence  upon  variables  to  be 
regulated  in  the  control  problem  and  require  substantially  fewer 
actuators  than  controlled  modes.  The  selection  process  itself,  as  well 
as  the  selections  it  produces,  exhibit  potentially  beneficial  design 
flexibility.  Linear-quadratic  disturbance-rejection  controllers 
employing  transducer  selections  generated  by  this  process  (cf.  Sec.  3) 
demonstrate  stability  in  the  presence  of  residual  modes  over  a  frequency 
range  substantially  beyond  that  of  the  design  modes.  These 


2-30 


demonstrations  are  significant  in  that  stability-enhancing  adjustments  of 
the  controller  feedback  structure  have  been  deliberately  postponed  in 
order  to  focus  on  the  connections  between  the,  selection  of  reduced-order 
models  and  active  transducers  on  the  one  hand  and  closed-loop  stability 
and  performance  on  the  other  hand. 

Subsequent  work  in  regard  to  active  transducer  selection  will 
include: 

(1)  Expansion  of  the  actuator  candidate  class  to  include 
rotational  actuators; 

(2)  Deeper  investigation  of  the  design  flexibility  inherent  in 
the  selection  process;  and 

(3)  Examination  of  the  interaction  between  the  selection  of 
reduced-order  models  and  the  selection  of  active  transducers. 

2.5  Appendix 

Theorem  2-6  is  an  immediate  consequence  of  the  following  two 
results. 

Proposition  2-7  (Column  Rearrangements)..  The  minimum  value  defined  by 
Eq.  (2-14)  is  unchanged  by  an  arbitrary  column  permutation  P:  mxm  of  A; 
i.e. , 

min  1  __  .  .  min  ... 

J(P  x  ;  AP,  b)  =  J(x;  A,  b)  (2-50) 

x  x 


proof.  The  identity 


J(P-1x; 


AP,  b)  = 


Jb-(AP) (P-1x) a2 


J(x;  A,  b) 


(2-51) 
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holds  for  all  n-vectors  x.  It  follows  that  for  each  minimum,  say  x0 
(cf.  Eq.  (2-13)),  of  J( • ;  A,  b) ,  P_1Xq  is  a  minimum  of  J( •;  AP,  b) , 
and  conversely.  Hie  desired  coincidence  of  minimum  values  Eq.  (2-50) 
follows  from  Eq.  (2-51)  by  this  correspondence  of  minimum  points.  □ 

Proposition  2-8  (Terminal  Removals).  Let  p  be  an  integer,  1  £  p  <_  m-1 . 
The  minimum  value  defined  by  Eq.  (2-14)  does  r.ot  decrease  if  columns 
p+1,»»»,m  are  removed  from  A;  i.e., 


min 


min 

xeE 


m  J(x;  A,  b) 


,  min 


1 ..*p' 


b) 


where  Al4.,p  denotes  the  submatrix  consisting  of  the  first  p  columns 
of  A. 

Proof.  Integers  k,  m,  and  n  are  defined  as  in  Eq.  (2-1).  In  view  of 
Proposition  2-7,  there  is  no  loss  of  generality  in  assuming. that  the 
leading  k  columns  are  linearly  independent.  Two  cases  are  considered, 
depending  on  whether  Eq.  (2-1)  represents  an  equality  or  a  strict 
inequality. 

Case  A;  k  a  min{n,m}.  If  n  <_  m,  then  the  result  is  trivial:  Eq.  (2-15) 
is  true,  and  therefore 


J  .  =  0  <  Jp. 

min  —  min 


Thus,  consider  that  m  <  n.  The  QR-decomposition  of  A  (cf .  Theorem  2-1 ) 
has  the  structure 


where  in  this  case  the  column  permutation  P  is  taken  to  be  the  identity 
matrix.  The  truncated  matrix  A^  . .  ,p  inherits  a  QR-decomposition  from 
Eq.  (2-52)  by  a  corresponding  truncation  of  R: 


S...P  -  5  (2t 

where  the  matrix  Q  is  the  same  as  in  Eq.  (2-52)  but  is  repartitioned  to 
be  compatible  with  (Ri)i...p:  pxp.  Denoting  by  q3  the  j-th  column 
of  Q,  it  can  be  seen  that 


(R, ) 


Using  Eq.  (2-19)  with  both  Q2  and  Q^,  it  follows  that: 


A  *  2 
IQ2  bj 


y.k 

^=p+i 


j (qj )*b| 2  +  KQ2*b|2 


> 


J  . 
min 


Case  B:  k  <  min  {n,m}.  The  proof  for  this  case  uses  the  same  approach 
as  in  Case  A,  except  that  Eq.  (2-52)  is  replaced  by  the  decomposition  for 
the  rank-deficient  case  that  is  a  simple  extension  of  Theorem  2-1: 


A  =  QR 


The  details  are  omitted.  □ 
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Table  2-2.  Reduced-order  models  in  terms  of  open-loop 
undamped  characteristic  frequencies: 

O  Interlaced  design 
•  Contiguous  design 
■  Contiguous  evaluation 


EIGENVALUE 


REAL  EIGENVALUES 


0.0 

0.0 

0.0 

0.0 

0.0 

0.0 _ 

8.650072E-01 
3.148905E400 
4.009007E+00 
4.434997E400 
8.638482E400 
1.344033E40X 
1.42424XE40X 
1.790326E401 
3.639226E401 
4.700934E401 
1.334769E402 
I.34171SE402 
1.408959E402 
1.564057E402 
1.675513E402 
2.372882E402 
2.412209E402 
4.149673E402 
1.051814E403 
1.055404E403 
2.449482E403 
2.4749S2E403 
3.037863E403 
3.040462E403 
3.083413E403 
3.183590E403 
3.853264E403 
4.265980E403 
5.379633E403 
5.939781E403 
7.075340E403 
7.319027E403 
6.&X4004E4Q3 
1 .065900E+04 
1.713255E404 
1.674883E404 
1.883882E+04 
1.884647E404 
1.884669E404 
1.889893E404 
2.227858E404 
2.364327E404 
2.455703E404 
2.643350E404 


0.0 

0.0 

0.0 

0.0 

0.0 

0.0 _ 

9.300699E-01 

1.774515E400 

2.00225IE4Q0 

2.105943E400 

2.939130E400 

3.666106E+00 

3.773911E400 

4.231224E400 

6.032599E40C 

6.862167E400 

1.1S5322E401 

I.I58326E401 

1.186994E401 

1.250623E401 

X.294416E401 

1.S40416E40I 

I.553129E401 

2.037074E401 

3.243167E40X 

3.248690E40X 

4.949223E401 

4.974889E401 

5.51X682E+01 

5.514038E+01 

5.552849E401 

5.642339E+01 

8.044223E40X 

6.53X447E40J. 

7.334598E401 

7.706998E+01 

0.41X504E40X 

8.S55X32E401 

9.281166E40X 

1.032424E402 

1.3089X4E402 

1.369264E402 

1.372546E402 

1.372824E+02 

X.372832E402 

1.374734E402 

1.49260XE402 

1.537637E402 

X.567068E402 

X.62S838E402 


0.0 

0.0 

0.0 

0.0 

0.0 

0.0 _ 

1.480284E-OX 

2.824228E-01 

3.186681E-01 

3.351712E-0X 

4.677770E-0X 

5.634789E-01 

6.006367E-0X 

6.734203E-0X 

9.601181E-0X 

1.092148E400 

1.8387S2E400 

X.843S33E400 

1.889X59E400 

1.990926E400 

2.060127E400 

2.451648E400 

2.471881E400 

3.242X05E40Q 

5.161660E400 

5.170464E400 

7.876933E400 

7.917781E400 

6.772115E400 

8.775864E400 

8.637633E400 

8.930061E400 

9.6X9679e»00 

1.0395X2E*OX 

1.167338E*0L 

1.226607E401 

1.338732E*0X 

1.36X592E401 

1.477143E401 

1.643153E401 

2.083200E40X 

2.179250E401 

2.134474E401 

2.184917E40X 

2.184930E401 

2.187956E401 

2.375S98E401 

2.447224E«0X 

2.494066E401 

2.58760XE401 
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LOSX:  61 
LOSY.  3E-5 


Figure  2-3.  Minimal  selection,  example  C. 
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LOSX:  0 
LOSY:  0 


Figure  2-5.  Preselection  D-17,  example  D. 


Figure  2-6.  Preselection  D-14,  example  D. 


m  »  1 


Figure  2-7.  Preselection  D-ll,  example  D 


Figure  2-8.  Minimal  selection  D-9A,  example  D. 


SECTION  3 


PROGRESS  ON  SYNTHESIS  OF  ACTIVE  CONTROL  FOR  BROADBAND 
DISTURBANCE  ACCOMMODATION 

PART  2s  CONTROL  SYSTEM  DESIGN  AND  EVALUATION 


3.1  Introduction 

The  philosophy  of  our  scientific  experiment  for  active  control 
design  with  ACOSS  Model  No.  2  was  reviewed  in  the  previous  section.  Now 
we  make  use  of  the  18  mode  (modes  7-24)  design  model  from  Revision  1  of 
Model  No.  2,  and  let  the  several  actuator  and  sensor  placements  he 
variables  in  the  experiment. 

This  section  will  review  the  control  system  design  philosophy 
which  remains  fixed  while  working  with  the  various  actuator  and  sensor 
selections.  After  the  design  results  are  examined,  some  changes  in  the 
controller  design  will  be  suggested  for  the  future  part  of  the 
experiment. 

3.2  Control  System  Design  Philosophy 

For  these  designs  we  have  employed  the  Linear-Quadratic-Gaussian 
(LQG)  methodology,  where  a  model  of  the  disturbance  has  been  augmented 
with  the  state  equations  of  the  structure.  Such  a  disturbance  model  uses 
knowledge  of  the  point  of  application  on  the  structure  and  power  spectral 
density  of  the  incoming  disturbances.  The  resulting  controller  is 
frequently  referred  to  as  a  disturbance  rejection  controller,  and  has 
been  treated  in  References  3-1  (Section  4)  and  3-2. 
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To  clarify  some  of  the  values  used  here  and  in  the  designs  of 
Reference  3-1,  we  examine  the  transfer  function  model  for  the  power 
spectral  density  function  corresponding  to  the  incoming  disturbances  as 
given  by  Eg.  (2-1)  of  Reference  3-1  (Section  2): 
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The  corresponding  state  space  representation  is  (cf.  Eq.  (2-28) ) 

v  =  BqUD 

where  up  is  the  solution  of 


u„  =  A  u„  +  B  oj 
D  L)  D  0) 

In  particular,  v  is  the  disturbance  acting  on  the  structure,  and  BD, 
the  influence  matrix  of  the  disturbances,  has  been  discussed  in  Section  2 
of  Reference  3-1.  Here 
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where  Xi  =  X2  ®  -u^  =  -31.4159,  and 
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where  =  /<5  uic  =  198.692.  The  above  terms  are  discussed 

in  Section  4  of  Reference  3-1,  but  without  noting  the  numerical  values. 
Also  the  covariance  matrix  of  the  fictitious  white  noise  input  vector  u 
is 


Cov(o)  oiT)  =  *2x2 


Note  that  the  power  spectral  density  function  has  a  bandwidth  of  5 
Hz  and  has  a  zero  frequency  value  of  40  N2/Hz  if  it  is  considered  to  be 
a  two-sided  function  of  frequency,  and  has  a  zero  frequency  value  of  80 
N2/Hz  if  it  is  considered  to  be  a  one-sided  function  of  frequency. 

Among  the  controller  designs  presented  in  this  section  for  the 
four  actuator/sensor  placements,  the  only  parameters  that  have  been 
altered  are  the  diagonal  terms  of  the  3x3  Q  matrix  in  the  performance 
index  for  the  LQ  regulator  where. 
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and  q1  weights  the  line-of-sigh t  displacements  in  the  x  direction 
(LOSX),  q2  weights  the  line-of- sight  displacements  in  the  y  direction 
(LOSY),  and  qj  weights  the  defocus  displacements  (Z).  In  particular, 
the  diagonal  terms  of  the  Q  matrix  are  adjusted  to  ensure  that  the  stare 
feedback  portion  of  each  LQG  design  will  meet  the  line-of-sight  (LOS) 
specifications.  These  LOS  specifications  are  taken  to  be: 
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Defocus  =  Z  =  [  £  (ui  -  Y3)j  /  _<  0.025  x  10  meters 
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where  Ui-Yj  represents  the  mean  square  ^alue  at  one  of  the  three  line- 
of-sight  outputs  resulting  from  the  i-th  disturbance.  For  our  designs 
there  are  two  identical  disturbances/  one  at  node  37  and  one  at  node  46 
of  the  model,  so  m  is  two  in  the  above  LOS  specifications. 

For  all  designs  the  control  weighting  matrix,  R,  in  the  LQ 
performance  index  was  an  identity  matrix  of  appropriate  dimension  and  its 
scalar  factor,  pR,  was  always  1x10"17. 

The  design  of  the  optimal  state  estimator  is  dual  to  that  of  the 
LQ  regulator  and  the  covariance  of  the  sensor  noise  term,  9,  was  taken  to 
be: 


Cov(8  9  )  =  pi 


o  mxm 


-17 

where  pQ  =  1x10  ,  and  m  is  the  number  of  sensor  outputs,  as  in  the 

designs  of  Section  4  of  Reference  3-1. 


The  state  estimator  is  connected  to  the  constant  LQ  gain  to  form  a 
compensator  which,  when  evaluated  with  the  18  mode  design  model  of  the 
structure,  is  guaranteed  to  yield  a  stable  closed-loop  system.  Note  that 
Model  No.  2  has  132  elastic  modes  of  higher  frequency  than  the  18  modes 
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used  in  our  design.  One  goal  should  be  to  have  a  control  system  robust 
enough  to  remain  stable  in  the  presence  of  these  132  high  frequency 
modes,  realizing  that  degradation  in  LOS  performance  would  occur.  Our 
designs  have  not  achieved  that  goal  at  this  point  in  the  experiment,  but 
the  results,  discussed  in  the  next  section,  are  encouraging. 

3.3  Controller  Results 

The  results  of  the  control  system  designs  with  four  of  the 
actuator/sensor  selections  described  in  Section  2  are  summarized  in  Table 
3-1.  The  diagonal  elements  of  the  Q  matrices  used  in  the  LQ  designs,  and 
the  resulting  line-of-sight  results  for  the  full-state  feedback  case 
under  the  column  headed  by  LQ  are  shown.  Note  that  the  LQ  designs  with 
translational  actuators/sensors  required  much  lower  weighting  of  the 
controlled  variables  (line-of-sight  displacements)  in  the  performance 
index  than  did  the  design  with  axial  members,  while  achieving 
approximately  the  same  LOS.  performance. 

Not  reflected  in  the  table  are  the  design  iterations  on  the  Q 
matrix  for  each  particular  actuator /sensor  selection  that  led  to 
acceptable  LOS  performance.  These  iterations  showed  that  each  Z  result 
was  relatively  insensitive  to  changes  in  its  own  weighting  or  the 
weighting  of  the  other  terms.  However,  an  order  of  magnitude  change  in 
the  weighting  for  LOSX  or  LOSY,  namely  qi  or  q2»  from  those  shown  in 
the  table  yielded  approximately  an  order  of  magnitude  change  in  the  LOSX 
or  LOSY  value. 

Recalling  that  actuator/sensor  selection  D-9A  was  a  subset  of 
selection  D-14,  one  can  note  only  a  slight  degradation  in  the  LOS 
responses  between  these  two  for  the  LQ  design. 

The  original  18  mode  LQ  design  for  selection  D-9A  encountered 
numerical  difficulties  when  evaluating  the  poles  of  its  s tate- feedback 
design.  These  numerical  difficulties  were  alleviated  when  mode  19,  which 
has  extremely  low  control  and  output  influences,  was  deleted  from  the 
design  process.  Mode  19  was  of  course  included  in  the  evaluation  models. 
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Going  to  the  final  column  of  Table  3-1  ,  one  can  see  that 
actuator/sensor  selection  C  yielded  the  largest  number  of  modes  in  a 
stable  evaluation.  The  addition  of  one  mode  beyond  any  of  the  stable 
evaluations  shown  in  the  last  column  gave  an  unstable  system.  Note  that 
selection  D-9A,  which  was  a  subset  of  selection  D-14,  had  a  much  lower 
stable  evaluation  model  than  D-14.  In  general,  taking  into  account  the 
number  of  actuators/sensors  in  each  selection,  the  translational  members 
yielded  control  systems  with  superior  stability  properties  compared  to 
the  design  using  axial  members. 

For  the  two  designs  giving  the  greatest  envelope  of  stability,  the 
LOS  results  were  computed  for  the  LQG  compensator  evaluated  with  the  18 
mode  design  model  and  with  a  30  mode  evaluation  model  (modes  7-36). 

While  slight  degradation  occurred  in  the  LCSX.  and  LOSY  results  for  the 
LQG  compensator  evaluated  with  18  modes  compared  to  the  full-state 
feedback  results  in  the  LQ  column,  significant  degradation  in  LOS  did 
occur  when  the  LQG  compensator  was  evaluated  with  modes  7-36. 

♦ 

3.4  Future  Controller  Designs 

As  part  of  the  scientific  experiment  explained  in  Section  2, 
changes  will  be  made  in  the  control  system  design  with  an  eye  toward 
improving  the  stability  results.  Specifically,  the  choice  of  state  and 
control  weighting  matrices  in  the  LQ  performance  index  could  be  altered 
to  lower  the  multivariable  bandwidth.  In  addition,  other  frequency 
shaping  techniques  may  have  to  be  used  to  ensure  adequate  high  frequency 
attenuation  in  the  LQ  design.  Hopefully,  in  addition,  a  reasonable 
disturbance  rejection  controller  design  can  be  found  which  relies  less  on 
explicit  knowledge  of  the  disturbances'  structural  locations  and 
stochastic  properties.  If  a  linear  quadratic  regulator  can  be  designed 
with  some  of  the  improved  properties  mentioned  above,  the  estimator  may 
have  to  be  designed  using  "full-state  recovery  techniques"  (Ref.  3-3)  to 
preserve  those  good  features,  which  otherwise  are  not  guaranteed  to  hold 
when  an  estimator  of  general  design  is  used  in  the  compensator. 
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Conclusion 


This  section  has  illustrated  the  LQG  disturbance  rejection 
controller  results  for  four  actuator  and  sensor  selections  used  on  ACOSS 
Model  No.  2.  These  selections  were  the  variables  under  investigation  in 
the  experiment  referred  to  in  Section  2.  The  results  indicated  that 
translational  actuators  and  sensors  gave  superior  line-of -sight  and 
stability  results  when  compared  to  axial  actuators  and  sensors.  Changes 
in  the  controller  design  philosophy  will  be  made  to  improve  the  stability 
results  in  the  future  part  of  this  experiment. 
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SECTION  4 


SAMPLED  DATA  CONTROL  OF  FLEXIBLE  STRUCTURES 
USING  NON-COLOCATED  VELOCITY  FEEDBACK 


4.1  Overview 

A  framework  is  developed  for  sampled  data  control  of  flexible 
structures,  in  terms  of  discrete  time  recursive  equations  in  second  order 
form.  This  framework  is  used  to  analyze  the  sampled  data  control  scheme 
where  the  loop  is  closed  using  constant  gain  output  velocity  feedback. 

It  is  well  known  that  the  closed  loop  is  stable  if  colocated  velocity 
feedback  with  symmetric  and  positive  definite  feedback  gain  is  used,  so 
long  as  the  sampling  rate  is  sufficiently  high.  In  this  section  it  is 
shown  that  the  closed  loop  can  be  stabilized  using  sampled  data  output 
velocity  feedback  for  arbitrary  sampling  rate;  our  approach  is  to  make 
use  of  appropriately  defined  non-colocated  velocity  feedback.  This 
approach  leads  to  explicit  stability  conditions  in  terms  of  the  feedback 
gain  and  the  sampling  rate. 

4 • 2  Models  for  Sampled  Data  Controlled  Flexible  Structures 

A  sampled  data  controlled  flexible  structure  can  be  defined  as  a 
distributed  parameter  system,  where  the  structure  input  is  the  output  of 
an  ideal  zero  order  hold  and  the  structure  output  is  sampled.  Although 
distributed  parameter  models  typically  involve  infinite  dimensional 
variables,  our  analysis  is  based  on  the  finite  dimensional  model 

Mq  +  Kq  =  Bu  (4-1  ) 
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For  simplicity  in  the  subsequent  development  no  structural  damping  is 
included.  The  structural  displacement  vector  q  =  (qi,»»»,qn)  and  the 
force  input  vector  u  =  (ui ,  •  •  *,1^) .  The  mass  matrix  M  and  the 
structural  stiffness  matrix  K  are  assumed  symmetric  and  positive 
definite.  Throughout,  we  consider  velocity  output  of  the  form 

y  =  Cq  (4-2) 

where  the  velocity  output  vector  y  =  (yi/**»/ym).  The  input 
influence  matrix  B  and  output  influence  matrix  C  are  assumed  to  be 
dimensionless. 

The  structure  input  u  is  defined  in  terms  of  the  input  sequence 
u^  by  the  ideal  zero  order  hold  relation 

u(t)  =  u  ,  kT  t  <  kT+T  (4-3) 

K  1 

The  output  sequence  y^  is  defined  in  terms  of  the  structure,  output  y  by 
the  ideal  sampling  relation 


yR  =  y(kT) 

The  fixed  value  T  >  0  is  the  constant  sampling  time.  This  open  loop 
sampled  data  controlled  structure  can  be  viewed  as  a  discrete  time  system 
with  input  sequence  and  output  sequence  y^,  where  k  =  0,1,* 

Let  $  be  a  nxn  nonsingular  modal  matrix  and  let  f)2  be  a  nxn 
diagonal  modal  frequency  matrix  [Ref.  4-1]  satisfying 

$TM$  =  I,  $'rK$  =  n2  =  diag(u)^,  •••,u/j) 


Introduce  the  coordinate  change 
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q 


=  $n 

so  that  Eqs.  (4-1),  (4-2)  can  be  written  as 

n  +  =  *Tbu  (4-4) 

y  =  C4>  n 

It  is  an  easy  task  to  solve  the  vector  equation  Eq.  (4-4),  using  Eq, 
(4-3),  to  obtain 

\+1  =  (cosOT)\  +  if1  (sinffT)^  +  ft“2(I  -  cosQT)  ^Bu^  (4-5) 

•  .-IT 

nk+1  =  -nfsinflTJr^  +  (cosffDr^  +  8  (sinflT)$  Bu^  (4-6) 

where  n  »  n(kT),  nk  =  n < kT ) ,  and  sinftT  =  diag(sin(<d^ T) ,  •  •  *,sin( u  T) ) , 
cosOT  =  diag(cos((UiT) •,cos({onT) )  (Ref.  4-2]. 

Although  the  first  order  recursive  equations  Eqs.  (4-5),  (4-6) 
could  be  used,  it  is  more  convenient  for  our  purposes  to  make  use  of  a 
second  order  recursive  equation  in  nk  alone 

nk+1  -  2 (cosflT) =  n"1  (sinflT)4>TB(uk  -  uk_1  )  (4-7) 

yk  =  C$  ^  (4-8) 
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The  modal  equations  Eqs.  (4-.7),  (4-8)  form  the  basis  for  our  subsequent 
analysis.  It  is  natural  to  make  use  of  the  recursive  equation  for 
alone  in  considering  velocity  feedback  systems. 

It  should  be  noted  that  relations  Eqs.  (4-7) ,  (4-8)  involve  no 
numerical  approximation;  they  are  valid  for  any  sampling  time  T  >  0. 

4.3  Constant  Gain  Velocity  Feedback 

Constant  gain  output  velocity  feedback  has  been  studied 
extensively  for  analog  controlled  structures.  Our  interest  is  in  use  of 
constant  gain  output  velocity  feedback  for  sampled  data  controlled 
structures. 

Consider  the  closed  loop  sampled  data  controlled  structure  defined 
by  Eqs.  (4-7)  and  (4-8),  using  the  control  input  sequence 

Uk  =  "Gyk  (4"9) 

where  G  is  a  constant  mxm  feedback  gain  matrix.  Substituting  Eq.  (4-9) 
into  'Eqs.  (4-7)  and(4-8),  a  closed  loop  recursive  equation  is  obtained 

nfc+1  -  [2cos(lT  -  jf1  (sinOT)$TBGC<t]nk 

+  [l  -  n'1  lsinflT)*TBGC*]J»k-1  =  0  (4-10) 

The  closed  loop  characteristic  eq  •  ,on  can  be  written  as 

d(T,z)  =  1  -  (2cosflT  -  ft  1  (sinfiT)$TBGC$)z 

'<1  -  n"1  (sintfmTBGC$)  ]  = 


0 


The  objective  of  constant  gain  velocity  feedback  control  is  to  make  the 
closed  loop  as  described  by  Eq.  (4-10)  geometrically  stable,  i.e. ,  to 
make  the  closed  loop  characteristic  zeros  lie  inside  the  unit  disk. 

We  use  Eq.  (4-10)  as  the  basis  for  our  subsequent  analysis  of  the 
closed  loop.  If  sinftT  is  nonsingular,  the  following  implications  hold: 
if  r^+0  as  k+®,  then  necessarily  u^+0  as  k->®  and  r^+0  as  k  •*■<*>; 

consequently  x^+0  and  x^+0  as  k+». 


4.4  Colocated  Velocity  Feedback 

Recall  the  following  results  for  constant  gain  output  velocity 
analog  feedback  control  where  u  =  -Gy.  If  colocated  force  actuators  and 
velocity  sensors  are  selected  so  that  C  =  BT,  then  the  closed  loop 
analog  controlled  structure  is  asymptotically  stable  if  G  is  any 
symmetric,  positive  definite  matrix,  and  if  a  certain  controllability 
assumption  is  satisfied  (Refs.  4-3,  4-4,  4-5] .  Moreover,  this  result 
does  not  depend  on  the  particular  values  of  the  modal  frequencies  and 
modal  functions. 

We  first  mention  a  rather  obvious  result  that  if  the  sampled  d.^ta 
feedback  control  is  chosen  according  to  the  analog  feedback  theory,  then 
the  closed  loop  is  stable  for  sufficiently  small  sampling  time.  The 
brief  proof  is  included  for  completeness?  it  also  serves  as  an 
introduction  to  our  subsequent  development. 

Theorem  4-1 .  Assume  that 
.  (a)  C  =  Bt, 

(b)  the  matrix  pair  (f)2,  ^Tg)  j_s  completely  controllable,  and 

(c)  G  is  symmetric  and  positive  definite. 

Then  the  closed  loop  equation  Eq.  (4-10)  is  geometrically  stable  for 
sufficiently  small  sampling  time  T  >  0. 


Proof.  Consider  the  associated  polynomial 


p(T,  w) 


det  [2(1  +  cosOT  -  n 


-1 


+  2ft"1 (sinftT) $TBGC$ 


(sinOT)  <!>TBGC$)^-r— — 
4 

Tw  1 

— 2  +  2(1  -  cosftT) ] 


3y  assmption  the  polynomial  defined  by 


(4-1  1  ) 


lim  p(T,w)  1— 

T+0  T2 

has  all  zeros  in  the  left-half -plane;  hence  there' is  a  number  T  >  0  such 
that  p(T,w)  ha^  all  its  zeros  in  the  left-half -plane  for  0  <  T  <  T. 

Using  the  bilinear  transformation 


z 


1  + 


1  - 


Tw 

2 

Tw 
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(4-1  2) 


it  follows  that  the  zeros  of  d(T,z)  are  necessarily  inside  the  disk; 
hence  Eq.  (4-10)  is  stable.  □ 

This  result  has  limited  application  since  there  is  no  indication 
of  the  range  of  values  of  the  sampling  times,  relative  to  the  feedback 
gain  matrix,  required  for  closed  loop  stability.  In  Reference  4-6, 
conditions  are  developed  which,  in  principle,  characterize  a  range  of 
values  of  the  sampling  time  for  which  the  closed  loop  system  is  stable. 
Unfortunately,  the  conditions  depend  on  an  a  priori  computable  bound  on 
an  exponential  matrix;  computation  of  such  a  bound,  in  analytical  terms, 
is  not  considered  in  Reference  4-6.  Of  course  one  could  perform  a 
numerical  search,  based  on  the  characteristic  polynomial  d(T, z)  (or 
equivalently  p(T,w))  for  a  specific  case,  to  determine  a  range  of  values 
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of  the  sampling  time  for  which  the  closed  loop  is  stable.  However,  for 
the  case  of  colocated  velocity  feedback  there  are  no  known  explicit 
conditions,  in  terms  of  the  sampling  time  and  feedback  gain  matrix,  which 
guarantee  stability  of  the  closed  loop  sampled  data  system. 

4.5  Non-colocated  Velocity  Feedback 

We  now  present  the  main  result  of  the  paper:  a  set  of  explicit 
conditions  on  the  sampling  time  and  feedback  gain  matrix  for  which  the 
closed  loop  sampled  data  controlled  structure  is  stable.  The  key  idea  is 
to  suitably  modify  the  assumption  of  colocated  force  actuators  and 
velocity  sensors. 

Theorem  4-2.  Assume  that  sinQT  is  nonsingular  and 

(a)  C  =  BT<&(sinfiT)$f V1  (f1  , 

(b)  the  matrix  pair  .(  [I  +  cosOT]  “U  I  -  cosfJT] ,  <t>TB )  is 
completely  controllable, 

(c)  G  is  symmetric  and  positive  definite,  and 

(d)  I  +  cosflT  -  T<tTCTGC4  is  positive  definite. 

Then  the  closed  loop  equation  Eq.  (4-10)  is  geometrically  stable. 

Proof.  The  assumptions  can  be  shown  to  guarantee  that  the  zeros  of 
p(T,w)  defined  in  Eq.  (4-11)  are  in  the  lef t-half-plane.  The  bilinear 
transformation  defined  in  Eq.  (4-12)  guarantees  that  the  zeros  of  d(T,z) 
are  necessarily  inside  the  unit  disk.  Hence  Eq.  (4-10)  is  stable.  □ 

Informally,  note  that  as  the  sampling  time  satisfies  T-*0  that  from 
condition  (a),  C-*-BT,  and  that  conditions  (b),  (c),  (d)  of  Theorem  4-2 
tend  toward  conditions  (b),  (c)  of  Theorem  4-1.  But  for  T  >  0  the  extra 
requirement  of  condition  (d)  of  Theorem  4-2  is  imposed  on  the  feedback 
gain;  this  imposes  a  limit  on  the  "size"  of  the  gain  matrix.  Condition 
(a)  of  Theorem  4-2  implies  that  the  output  velocity  feedback  given  by 
Eq.  (4-2)  be  selected  in  a  specific  way.  Clearly  the  actuators  and 
sensors  would  generally  not  be  colocated. 
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Conditions  (a)  and  (d)  of  Theorem  4-2  do  depend  on  explicit 
knowledge  of  the  modal  frequencies  and  modal  functions  and  they  depend  on 
the  particular  value  of  the  sampling  time.  Unlike  the  corresponding 
analog  case,  the  closed  loop,  with  assumptions  of  'Theorem  4-2  satisfied, 
is  not  robust  to  arbitrary  uncertainties  in  the  modal  data  and  to  modal 
truncation.  However,  an  arbitrary  structure  can  be  stabilized  using 
sampled  data  velocity  feedback  with  arbitrary  sampling  time;  the 
suggested  approach  depends  on  a  proper  choice  of  force  actuator  and 
velocity  sensor  locations  so  that  condition  (a)  of  Theorem  4-2  is 
satisfied. 

4.6  An  Example  • 

As  a  simple  example  consider  the  scalar  equations  in  modal 
coordinate  form  with  one  actuator  input 

"  2 

n1  +  a,  n,  =  u 

n2  +  4*2  =  u 


2  2 

where 

Based  on  Theorem  4-1 ,  the  sampled  data  feedback  control 


u ( t)  =  -G^fkT)  +  n2< kT)  ]  ,  kT  £  t  <  kT+T, 

with  G  >  0  and  w^T  *  nir,  n=0,1,««»;  i=1,2  stabilizes  the  closed  loop 
for  T  >  0  sufficiently  small.  Even  in  this  simple  case  analytical 
conditions  on  the  range  of  values  of  the  sampling  time  and  the  feedback 
gain  for  wl.ich  the  closed  loop  system  is  stable  are  exceedingly 
complicated.  Evidently  a  proper  choice  of  values  for  the  feedback  gain 
and  the  sampling  time  would  require  use  of  numerical  search  procedures. 
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Now  consider  the  sampled  data  control 


sinco.T  sinio-T 

u(t)  3  -sfC  -^r~  )  VkT)  +  )  ^2(kT)]'  ItT  <.  t  <  kT+T, 

where  the  output  is  chosen  so  that  condition  (a)  of  Theorem  4-2  is 
satisfied;  clearly  such  velocity  sensors  would  not  be  colocated  with  the 
specified  force  actuators.  On  the  basis  of  Theorem  4-2  the  closed  loop 
is  guaranteed  to  be  stable  if  ojj,T  t  nir,  n=0,1,«*«;  i=1,2;  and  if  the 
relatively  simple  analytical  conditions 

2 

(ui.T)  (1  +  cosoj.T) 

0  <  GT  <  — - - — - —  ,  i=1,2 

(sinw^T) 

2  2 

„  {(i).  T)  ( U)  T)  {1  +  cosoj.T)  (1  +  cosoj_T) 

U  <  CaT  \ _ ] _ Z _ 1 _ Z _ _ 

2  ’22  2 
(oj2T)  (1  +  cos^T)  (sim^T)  +  (o^T)  (i  +  cosoj^T)  (sinu)2T) 

are  satisfied.  Thus  proper  choice  of  values  for  the  feedback  gain  and 
the  sampling  time  is  considerably  simplified. 

4.7  Conclusions 

We  have  presented  one  approach  where  feedback  gains  and  sampling 
time  can  be  readily  selected  to  guarantee  that  the  sampled  data 
controlled  structure  is  stable.  Our  conditions  imply  that  force 
actuators  and  velocity  sensors  be  located  to  satisfy  a  specific 
analytical  condition  which,  in  fact,  depends  on  the  sampling  time.  These 
explicit  results  are  in  contrast  to  the  lack  of  available  explicit 
guidelines  for  choosing  the  feedback  gains  and  a  sampling  time  in  the 
case  where  colocated  actuators  and  sensors  are  used. 
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SECTION  5 


OPTIMAL  TRACKING  AND  TERMINAL  TRACKING  MANEUVERS 
FOR  FLEXIBLE  SPACECRAFT 


5.1  Introduc tion 

The  problem  of  feedback  control  of  flexible  spacecraft  undergoing 
large-angle  maneuvers  is  a  topic  of  continuing  interest  in  the  aerospace 
community.  Two  types  of  closed-loop  slewing  maneuvers  are  of  particular 
interest.  First,  tracking  maneuvers  where  the  control  is  determined  in 
such  a  way  as  to  cause  the  state  to  track  or  follow  a  desired  output 
state.  Second,  terminal  trackin'.,’  maneuvers  where  some  elements  of  the 
state  are  required  to  exactly  satisfy  nonzero  terminal  constraints. 

The  necessary  conditions  defining  the  solutions  for  tracking  and 
terminal  tracking  problems  are  well  known.  Nevertheless,  the 
computational  burden  required  to  numerically  integrate  the  resulting 
systems  of  coupled  nonlinear  Riccati-like  differential  equations  has 
seriously  hindered  the  common  application  of  these  control  techniques. 

To  overcome  the  computational  difficulties  noted  above,  closed 
form  solutions  are  presented  for  the  necessary  conditions  defining  the 
solutions  for  the  following  linear  time-invariant  control  problems: 

(a)  The  linear  tracking  problem; 

(b)  The  terminal  tracking  problem; 

(c)  The  state  trajectories  for  the  conventional  closed-loop 
system  dynamics  equation. 

The  remainder  of  the  section  is  presented  in  three  parts.  Section 

5.2  presents  closed  form  solutions  for  the  Riccati  and  prefilter 
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equations  of  the  linear  tracking  problem.  Maneuver  simulations  are 
presented  for  both  conventional  and  control-rate-penalty  control  designs 
[Refs.  5-1,  5-2].  Section  5.3  presents  closed  form  solutions  for  the 
three  nonlinear  coupled  Riccati-like  differential  equations,  defining  the 
solution  for  the  terminal  control  problem.  Maneuver  simulations  are 
presented  for  both  conventional  and  control-rate-penalty  control 
designs-.  Section  5.4  presents  a  closed  form  solution  for  the  state 
trajectories  of  a  conventional  feedback  control  system,  when  the  plant 
and  the  state  are  assumed  to  be  perfectly  known.  The  solution  for  the 
state  is  presented  in  a  simple  recursive  form. 


5.2  Closed  Form  Solution  for  the  Linear  Tracking  Problem 

In  this  subsection  the  problem  of  maneuvering  a  flexible 
spacecraft  through  a  large  angle  is  considered,  where  the  feedback 
control  system  is  required  to  track  a  desired  output  state.  The  desired 
output  state  is  assumed  to  be  provided  from  an  open-loop  solution  for  the 
linear  time  invariant  system  model  (see  References  5-1  and  5-2  for 
examples  of  suitable  open-loop  solutions). 

The  optimal  control  problem  of  this  section  is  formulated  by 
defining  a  performance  index  which  consists  of  an  integral  of  quadratic 
forms  in  the  state,  control,  and  control-rate.  The  necessary  conditions 
for  this  problem  lead  to  two  nonlinear  Riccati-like  differential 
equations  [Ref.  5-6].  One  equation  is  the  standard  Riccati  equation, 
which  possesses  a  well-known  solution  in  terms  of  a  steady-state  plus 
transient  term  (Ref.  5-1].  The  remaining  vector  equation  is  a  prefilter 
equation  which  has  foreknowledge  of  the  desired  output  state,  and  is 
coupled  to  the  standard  Riccati  equation.  In  Section  5.2.4,  a  new  closed 
form  solution  is  presented  for  the  auxiliary  Riccati-like  vector 
equation. 
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5.2.1  Optimal  Linear  Tracking  Problem 

The  optimal  reference  tracking  problem  is  formulated  by  finding 
the  control  inputs  u(t)  to  minimize 

•J  =  y  lz*(tf) 


Z(V«S+  2 


if'  [« 


z*(t)  -  z(t)J^  +  liu(t)  n^]dt 


(5-1) 


for  the  system 


x  =  Ax  +  Bu  ,  x(t  )  =  x  (5-2) 

o  o 

z  =  Cx  (5-3) 

T 

where  1  p |  =  p  Wp,  z*  is  the  desired  output  state,  A  is  the  system 

w 

dynamics  matrix,  B  is  the  control  influence  matrix,  C  is  tne  measurement 
influence  matrix,  S  is  the  terminal  state  weight  matrix,  Q  is  the  state 
weight  matrix,  R  is  the  control  weight  matrix,  x  is  the  state,  and  u  is 
the  control.  (For  maneuvers  where  the  state  is  augmented  by  the  control 
and  control-rate  penalties,  the  A  and  B  matrices  of  Eq.  (5-2)  are 
modified  as  shown  in  Reference  5-1). 

5. 2. 1.1  The  Open-Loop  Control  Problem  for  the  Desired  Output  State 

The  desired  output  state,  z*(t),  in  Eq.  (5-1)  is  assumed  to  be 
provided  by  the  solution  for  the  following  open-loop  control  problem, 
where  we  seek  the  control  inputs  u*(t)  to  minimize 

1  tf  ?  ? 

J  *  -■  J  J  +  Hu*  (t)  lRjk  ]dt  (5-4) 

*o 
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for  the  system 


x*  =  Ax*  +  Bu*,  x*(t  )  = 

o 


X*,  X*(tf)  =  X* 


(5-5) 


z*  =  Cx* 

As  shown  in  Reference  5-2,  the  closed  form  solution  for  the 
open-loop  problem  defined  by  Eqs.  (5-4)  and  (5-5)  can  be  written  as 


where 


i  ^  ^  | 

U*(t)/ 


Q(t-t  ) 
e  o 


x* 

o 


X* 

o , 


(5-6) 


A 

L  -Q* 


-1  T 

-B(R* )  B 


-A 


(5-7) 


Xo 


(*  y 

■’/x*  -  *  x*l 

\  x*X*  / 

l  ,x*x*  ; 

(5-8) 


( .) 

e  =  the  exponential  matrix, 

_X*(t)  =  the  costate  vector  for  the  open-loop 
solution, 


and  ^x*X*  are  Part^t^ons  f°r  t=tf* 


To  obtain  x*(t)  from  Eq.  (5-6),  the  following  mapping  equation  is 
defined: 
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* 


X*(t) 


X* 

o 


(5-9) 


=  H  e 


n(t- 


t  ) 
o 


X* 

o 


where  H  =  [i  o]  is  a  (n  x  2n)  selection  operator  and  x*  is  (n  x  i ) . 

As  shown  in  what  follows,  the  explicit  presence  of  e<  *)  in  Eq. 

(5-9)  permits  a  closed  form  solution  to  be  obtained  for  the  prefilter 
equation  of  the  linear  tracking  problem. 

5.2.2  Necessary  Conditions  for  the  Linear  Tracking  Problem 

The  necessary  conditions  defining  the  optimal  tracking  problem  for 
Eqs.  (5-1),  (5-2),  and  (5-3)  are  gi/en  by  the  following  Riccati-like 
differential  equations  [Ref.  (5-3),  pp.  100-102]: 

fi  a  -PA  -  AT?  +  PBR-1  BTP  -  CTQC  ;  P(tf)  *  CTSC  (5-10) 


l  a  -[a  -  BR-1BTp]Tr  -  CTQCX*  }  5(tf)  a  CTSCxMtf)  (5-11) 

where  the  optimal  control  is  given  by 

u  *  -R_,BT[px  -  (5-12) 

5.2.3  Closed  Form  Solution  for  the  Time-Varying  Riccati  Equation 
The  solution  for  Eq.  (^-10)  has  been  previously  obtained  in 
References  5-2  and  5-4,  and  can  be  written  as 


P(t)  =  P  +  0_1 (t) 
ss 


(5-13) 
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where  Pss  is  the  solution  for  the  algebraic  Riccati  equation 


-P  A  -  ATP  +  P  BR-1BTP  -  CTOC  =  0  (5-14) 

ss  ss  ss  ss 


and  the  closed  form  solution  for  0(t)  can  be  shown  to  be 


e<t)  =  0ss  +  aA(t'tf)  [0f  -  0sg]eA  (t"V  (5-15) 


-  -IT  T  -1 

where  A  =  A  -  BR  BP  ,  0  ^{C  SC-P  )  ,  and  0  satisfies  the  algebraic 

ss  f  ss  ss 

Lyapunov  equation 


A0  +  0  AT  =  .BR-1BT  (5-16) 

SS  ss 

5.2.4  Closed  Form  Solution  for  the  Prefilter  Equation 

The  solution  for  5  in  Eq.  (5-11),  follows  on  assuming  the  product 
form  solution 


5(t)  =  Q~\t)  r(t)  (5-17) 

where  0  is  defined  by  Eq.  (5-15)  and  the  linear  constant  coefficient 
vector  differential  equation  for  r  can  be  shown  to  be 


r  -  Ar  =  -0CTQCx*  ,  r(tf)  =  0fCTSCx* (tf )  (5-18) 

(To  derive  Eq.  '.5-18)  one  requires  the  differential  equation  for  0  given 
bys  0  =  A0  +  ©AT  -  BR-1BT). 

The  solution  for  c  in  Eq.  (5-18;  can  oe  shown  to  be 

r(t)  »  eAt[rQ  -  /  e“AT0( t)CTQCx* ( r 1  dx]  (5-19) 

o 
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where  the  initial  condition  for  the  vector  r  is  given  by 


r 

o 


-At  t  — 

e  E  0CSCx*(t  )  +  /  e  T0(t)CTQCx*(t)  dT 


(5-20) 


The  solution  for  the  integral  expression  in  Eq.  (5-19)  follows  on 
introducing  0(t)  from  Eq.  (5-15)  and  x*(t)  from  Eq.  (5-9)  into  Eq. 
(5-19),  leading  to 


/  f"AT  0(t)CTQCx*(t)  dt 
o 


t 

r  -At  {It  , 

J  e  dT 


where 


/ 


033cTsch 


D 


2 


(5-22) 


(5-23) 


d3  =  ctqch 


(5-24) 


As  shown  in  Reference  5-5,  the  integrals  involving  the  matrix 
exponentials  in  Eq.  (5-21)  can  be  easily  evaluated  by  defining  the 
constant  matrices 


r  -  i 

•  -O' 

> 

a 

-A*  D3 

C1  - 

'  c2  = 

1 

o 

I _ 

- , 

a 

o 

_ i 
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and  computing  the  matrix  exponentials 


F1  (t)  G^ ( t) 

"F3(t)  G2(t)’ 

V 

C2fc 

e  = 

/  e 

.  0  F2(t)- 

.  0  v*>_ 

where 

F1  (t)  =  eAt,  F2<  t)  =  P4(t)  -  eQt,  F3(t)-  =  e-^* 

G^t)  =  eAt  /  e"ATDienT  dT,  G^t)  =  J*  8*%  eflT  dT 

o  o 


As  a  result,  it  follows  that  the  integrals  in  Eq.  (5-21)  are  given  by 


/  e"ATDienTdT  =  Fj ( t)G^  ( t) 


(5-25) 


t  -T 


/  eA  TD3enTdT  =  F^(t)G2(t) 


(5-26) 


5*2. 4.1  Recursion  Relationsnips  for  the  Integrals  of  Exponential 
Matrices 

In  order  to  evaluate  Eqs.  (5-25)  and  (5-26)  at  discrete  time 
steps,  the  semi -group  properties  of  exponential  matrices  are  now 
exploited;  yielding  the  following  recursion  relationships  for  the  matrix 
partitions  of  eC>(n+')4t  and  .C2(nH)4t 
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F.  [(n+1 )  At]  =  F  (At)  Fj  (n At)  ,  F_.  (0)  =  I  ,  j  =  1,2, 3, 4 

G  [(n+1  )At]  =  F  (At)G1  (nAt)  +  G  ( At)F  (nAt)  ,  G^O)  =  0 

G2  [ (n+1 ) At]  =  F3(At)G2(nAt)  +  G2( At)F2(nAt)  ,  G2(0)  =  0 

from  which  the  integrals  of  Eq.  (5-25)  and  (5-26)  follow  as 

(n+1) At  -  _  _ 

/  e  =  F*  [ (n+1 ) At]  Gjtn+DAt]  (5-27) 

o 


(n+1) At  -T 
/  eA  TD.efiTdT 


F^[(n+1 ) At ]  G2[(n+1)At] 


(5-28) 


where  At  =  (tf  -  tQ)/m  and  m  is  the  total  number  of  discrete  time 
steps. 


5. 2. 4. 2  Vector  Recursion  Relationships  for  r[(n+1)At] 

Substituting  Eq.  (5-27)  and  (5-28)  into  Eq.  (5-21),  yields 

(ntmt  -At  T 
/  e  0(t)C  QCxMtMt  3 

o 

=  F3[(n+1)At]  G1 [(n+1)At]s 

+  D2F^*[  (a+i  )  Ac]G0  [  (n+1 )  At  ]s 


or 
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F^[(n+1 )At]{F1 (At)G1 (nAt)  +  G1 ( At)F2(nAt) }sc 


+  D3F^[(n+1  )At]{F3(At)G2(nAt)  +  G2( At) F^nAt)  Js^ 


where  s 


[<*s>T  <^>T1 


The  number  of  mathematical  operations  required  in  the  equation 
above  can  be  minimized  by  defining  the  following  three  vector  functions 
and  associated  recursion  relationships: 

x'  f  (n+1 )  At  1  =  G.  [  (n+1 )  At  Is  =  F,  (  At)x,  (nAt)  +  G  ( At)x_  (nAt) ;  x  (0)  =0 

1  1  J  O  1  1  1  j  i 

x  f(n+1)At]  =  G  [ (n+1 ) At]s  =  F  (At)x  (nAt)  +  G  ( At)x3(nAt) ;  x2<0)  =  0 


[(n+1)At]  =  F2 [ (n+1 ) At]sQ  =  F2( At)x3 (nAt) ;  x3(0)  =  s 


As  a  result,  the  recursive  expression  for  the  integral  in  Eq. 
(5-19),  can  be  written  as 


(n+1 )  At 

/  e"nT0(T)C1QCx<'(T)dT 
o 


F3  [  (n+1 )  At]  x1  [  (n+1)  At] 


+  D3F^[(n+1 ) At ]  x2[(n+1 )At] 


where  the  recursion  relationships  for  the  matrix  partitions  G i [ ( n + 1 ) At ] 
and  G2[(n+1)At]  have  been  replaced  by  the  vector  recursion 
relationships  for  X}[(n+1)At],  X2[(n+1)At],  and  X3[(n+1)At]. 
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5. 2. 4. 3  Recursion  Relationship  for  Q[(n+1  )At] 


*  Subject  to  the  recursion  relationships  above,  the  time  histories 
for  P  and  £  in  Eq.  (5-12)  easily  follow,  where  o[(n+1)At]  in  Eq»  (5-15) 
is  given  by 


©[  (n+1 )  At]  =  0  +  F  (  At)  [©(n At)  -  0  jF^lAt) 

O  O  J  O  w  i 

and 

0(0)  -  eM  +  »J(tf)[e£  -  8sl)Fj(tf) 

5.2.5  Illustrative  Examples 

The  specific  model  considered  in  this  section  (see  Figs.  5-1  and 
5-2)  consists  of  a  rigid  hub  with  four  identical  elastic  appendages 
attached  symmetrically  about  the  central  hub.  In  particular,  the 
following  idealizations  are  considered:  (i)  single-axis  maneuvers;  (ii) 
in-plane  motion;  (iii)  anti-symmetric  deformations;  (iv)  small  linear 
flexural  deformations;  (v)  only  the  linear  time-invariant  form  of  the 
equation  of  motion  is  considered;  and  (vi)  the  control  actuators  are 
modelled  as  concentrated  torque  generating  devices.  The  distributed 
control  system  for  the  vehicle  consists  of:  (i)  a  single  controller  in 
the  rigid  part  of  the  structure,  and  (ii)  a  single  controller  in  each 
elastic  appendage.  In  the  optimal  control  performance  index,  however, 
the  control  weighting  matrices  are  adjusted  in  order  to  have  the  rigid 
body  control  provide  the  primary  torque  for  maneuvering  the  vehicle, 
while  the  appendage  controllers  act  principally  as  vibrational 
suppressors. 

For  all  cases,  the  following  configuration  parameters  are  assumed 
the  moment  of  inertia  of  fte  undeformed  structure  about  the  spin  axis,  I 
is  6764  kg-m^;  the  mass/length  ratio  of  the  four  identical  elastic 
appendages,  p,  is  0.04096  kg/ra;  the  length  of  each  cantilevered 
appendage,  L,  is  35  m;  the  flexural  rigidity  of  each  cantilevered 
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appendage,  El,  is  1500  kg-m3/s2;  and  the  radius  of  the  rigid  hub,  r, 
is  1  m.  For  simplicity,  each  appendage  is  assumed  to  have  one  controller 
located  half-way  along  its  span.  In  the  integrations  over  the  mass  and 
stiffness  distributions,  the  radius  of  the  hub  is  not  neglected  in 
comparison  to  the  appendage  length.  In  the  structural  modeling 
equations,  the  following  comparison  functions  have  been  adopted  as 
"assumed  modes" 


<jMx-r)  =>  1  -  cos[pir(x-r)/L]  +  0 . 5 ( — 1  )P+1  [prr(x-r)/L]' 


(p  =  1 , 2, •  •  •  ,n)  (5-29) 


which  satisfy  the  geometric  and  physical  boundary  conditions. 


4>  -4'  =  =  d,*  "  =0  (5-30) 

Tp  x-r  Tp  x~r  Tp  x=r+L  Tp  x=r+L 


of  a  clamped  free  appendage,  where  (*)'  =  d(*)/dx.  In  addition,  full 
state  feedback  is  assumed  for  the  results  of  this  section. 

The  output  vector  is  assumed  to  be  given  by 


y  -  {yt  yj  y3  y4  y*}*  (5-3 d 

where 

y^  «  central  hub  angular  velocity 


y2  =  four  tines  the  angular  velocity  at  each  apperdage  controller 
location  relative  to  the  central  hub 

y^  a  hub  angular  position 


y^  =  four  times  the  tip  deflection  of  each  appendage 


y,  =  controls  and  control-rates 
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The  number  of  elements  in  is  given  by  the  number  of  controllers  on 
each  of  the  four  appendages.  The  output  sub-vector  y^  exists  when 
control-rates  are  penalized.  The  elements  of  y^  correspond  to  the 
cidditionally  augmented  states. 

Referring  to'Tables  5-1  a  and  5-2,  the  graphical  summaries  of  the 
states  and  controls  are  discussed  qualitatively  in  what  follows. 

Case  1  (Fig.  5-3)  presents  a  terminal  controller  example  maneuver 
in  which  two  flexible  modes,  are  controlled  in  addition  to  the  rigid  body 
rotation.  There  are  six  states  in  this  problem,  all  of  which  are 
specified  at  the  final  time.  The  time  histories  for  this  case  are 
virtually  identical  to  the  results  obtained  from  an  open-loop  controller 
using  the  same  weight  matrices  in  the  performance  index.  The  jump 
discontinuities  in  the  control  torques  at  the  initial  and  final  times  are 
characteristic  of  controllers  with  no  penalties  on  the  control-rates  in 
the  performance  index.  As  shown  in  Reference  5-1,  the  presence  of  these 
discontinuities  tends  to  excite  both  modeled  and  unmodeled  high  frequency 
structural  modes. 

Case  2  (Fig.  5-4)  presents  the  same  maneuver  as  in  Case  1,  except 
that  the  performance  index  now  includes  penalties  on  the  first  and  second 
time  derivatives  of  the  control  torques  [Ref.  5-2].  Including  these 
penalties  allows  the  control  torques  and  torque-rates  to  be  specified  at 
the  initial  and  final  times,  thereby  eliminating  the  terminal  jump 
discontinuities  in  the  control  profiles,  ur  this  particular  control 
problem,  the  states  are  augmented  to  include  the  control  and  first  time- 
derivative  of  the  control,  while  the  second  time -derivative  of  the 
control  becomes  the  commanded  input  to  the  augmented  system.  On 

comparing  with  the  results  of  Case  1,  one  finds  much  smoother  modal 

• 

amplitude  and  control  torque  time  histories,  although  the  peak  torques 
and  tip  deflection  are  somewhat  higher  than  those  in  Case  1 . 
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Case  3  (Fig.  5-5)  presents  the  results  of  a  spin-to-rest  maneuver 
in  which  the  weight  matrices  and  final  conditions  are  identical  with 
those  of  Case  2.  The  initial  conditions  are  different  from  those  in  Case 
2,  but  the  time-varying  feedback  gains  are  identical.  As  the  plots  in 
Figure  5-7  show,  the  terminal  controller  has  little  difficulty  in 
bringing  the  system  to  the  required  final  conditions. 

Case  4  (Fig.  5-6)  presents  the  results  of  a  spin-to-rest  maneuver 
where  off-nominal  plant  parameters  are  used.  In  particular,  the 
moment  of  inertias  for  the  hub  and  appendages  are  increased  by  1 0%  and 
the  appendage  mass  per  unit  length  is  decreased  by  10%,  while  the 
time-varying  feedback  gains  are  computed  using  the  nominal  structural 
parameters.  Although  the  structural  parameters  are  changed  by  only  10%, 
the  overall  effect,  on  the  entire  structure  is  quite  large.  Specifically, 
the  first  and  second  mode  eigenvalues,  u)2,  are  increased  by  16%  and  22% 
respectively,  while  the  ratio  of  the  hub  moment  of  inertia  to  appendage 
moment  of  inertia  is  increased  by  22%.  Nevertheless,  the  plots  show  that 
most  of  the  state  variables  reach  their  prescribed  terminal  boundary 
conditions. 

In  Cases  1  through  4  the  results  have  been  compared  with  the 
corresponding  open-loop  maneuvers,  and  have  been  found  to  be  identical 
within  plotting  accuracy,  as  expected.  Furthermore,  the  closed  form 
solutions  for  P(t)  and  £(t)  have  alsc  been  compared  with  the  backward 
integrated  solutions  of  Eqs.  (5-10)  and  (5-11)  for  verification. 


5.3  Closed  Form  Solution  for  the  Terminal  Tracking  Problem 

In  this  subsection  the  problem  of  maneuvering  a  flexible 
spacecraft  through  a  large  angle  is  considered,  wh*re  the  terminal  state 
is  subject  to  constraints.  The  basic  approach  differs  from  the  standard 
linear  quadratic  regulator  design,  in  that  a  quadratic  penalty  term  on 
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the  final  value  of  the  state  is  not  included  in  the  performance  index. 

As  a  result,  the  numerical  difficulties  associated  with  handling  large 
terminal  penalty  matrices  are  not  encountered. 

The  optimal  control  problem  of  this  section  is  specified  by 
defining  a  performance  index  which  consists  of  an  integral  of  quadratic 
forms  in  the  state,  control,  and  control-rates,  where  the  terminal  state 
is  subject  to  constraints.  In  particular,  the'  terminal  constraints 
replace  the  conventional  weighted  quadratic  form  in  the  terminal  states. 
The  necessary  conditions  for  this  problem  lead  to  three  nonlinear 
Riccati-like  differential  equations  [Ref.  5-6] .  One  equation  is  the 
standard  Riccati  equation,  which  possesses  a  well-known  solution  in  terms 
of  a  steady-state  plus  transient  term  (see  Eq.  (5-13)).  On  the  other 
hand,  the  remaining  two  equations  are  coupled  to  the  standard  Riccati 
equation.  In  Sections  5.3.3  and  5.3.4,  new  closed  form  solutions  are 
presented  for  the  two  auxiliary  Riccati-like  equations. 


5.3.1  Optimal  Terminal  Controller 

The  optimal  terminal  control  problem  is  formulated  by  finding  the 
control  inputs  u(t)  to  minimize 

tf 

J  =  J  f  (*VqPx  +  uTRu)  dt  (5-32) 

t 

o 

for  the  system 

x  =  Ax  +  Bu  ,  given  x(t  )  (5-33) 

o 


with  outputs 
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y 


Fx 


(5-34) 


where  x,  u,  A,  B,  Q,  and  R  are  defined  following  Eq.  (5-3),  F  is  the 
measurement  influence  matrix,  and  the  stace  is  subject  to  the  specified 
terminal  constraints 


xi(tf)  =  x^^  ,  i  =  ),•••, q  (q  <_  n)  (5-35) 

Of  particular  interest  is  the  fact  that  the  performance  index  of  Eq. 
(5-32)  does  not  contain  a  terminal  weight  matrix  which  penalizes  the 
final  values  of  the  state. 

5 . 3 . 1 . 1  Necessary  Conditions  for  the  Terminal  Tracking  Problem 

As  shown  in  Reference  5-6,  the  necessary  conditions  defining  the 
optimal  solution  are  given  by  the  following  coupled  Riccati-like  matrix 
differential  equations 

P  +  PA  +  ATP  -  PBR-1BTP  +  FTQF  =  0  ;  P(tf)  =>  0  (5-36) 

S+  (AT  -  PBR_1BT)S  =  0  ;  ST(tf)  =  (3i|>/3x)|t  (5-37) 

G  a  STBR_1 BTS  ;  G(tf)  =  0  (5-38) 


where 
given  by 


[x, ,  x_,  •••,  x  ]  and  the  optimal  continuous  feedback  law  is 
1  2  q 


U(t)  =  -C(t)x( t)  -  Dtt)^ 


(5-39) 


where 
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C  =  R  B  (P  -  SG  S  ) 

_  -1  _T_  -1 

D  =  R  B  SG 


(5-40) 

(5-41) 


The  control  vector  u  of  Eq.  (5-39)  takes  the  state  vector  from  x(tQ)  at 
time  t0  to  at  time  tf,  while  minimizing  the  cost  function  of  Eq. 

(5-32)  . 


In  the  next  section  closed  form  solutions  are  presented  for  Bgs. 
(5-36),  (5-37),  and  (5-38),  thus  reducing  the  solution  for  the  feedback 
gains  to  the  direct  computation  of  algebraic  equations  without  numerical 
integration. 


5.3.2  Closed  Form  Solution  for  the  Time-Varying  Riccati  Equation 


P(t)  «  P_+  0-\t) 

ss 


(5-42) 


where  Pgs  is  the  solution  for  the  algebraic  Riccati  equation 

-ATP  -  P  A  +  P  BR-1BTP  -  FTQF  a  0  (5-43) 

ss  ss  ss  ss  * 


and  the  solution  for  0(t)  follows  as 


0(t) 


ss 


A(t-tf ) 
e 


(© 


ss 


+  p 


r1 

SS 


A  (t-t  ) 
)e 


(5-44) 


where  is  the  exponential  matrix,  A  =  A 

for  0ss  is  defined  by  Eq.  (5-16). 


-1  T 

BR  B  P  , 
ss 


and  the  solution 


5.3.3  Closed  Form  Solution  for  S(t) 

The  new  solution  for  S(t)  in  Eq.  (5-37)  follows  on  assuning  the 
product  form  solution 
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S  ( t ) 


(5-45) 


=  o-1 (t)S  (t) 

c 


where  0  is  defined  by  Eq.  (5-44)  and  Sc  is  to  be  determined. 

Substitution  of  Eq.  (5-45)-  into  Eq.  (5-37)  leads  to  the  following 
linear  constant  coefficient  matrix  differential  equation  for  Sc(t): 

S  (t)  -AS  (t)  «  Oj  S  (t_)  =  -P_1  S(t.)  (5-46) 

c  c  c  f  ss  r 

from  which  it  follows  that  the  solution  for  S(t)  in  Eq.  (5-45)  can  be 
written  as: 


A(t-tjf)  _■) 

S(t)  »  -0  ( t)e  P  S(t_) 

ss  t 

.  -T  -IT 

where  the  differential  equation  for  0  is  given  by:  0  =  A0  +  0A  •  -  BR  B  . 
The  solution  for  Sc(t)  in  Eq.  (5-46)  can  be  shown  to  be 

Mt-t  ) 

S  (t)  *  -e  P  S(t„)  (5-47) 

c  ss  t 

5.3.4  Closed  Form  Solution  for  G(t) 

The  solution  for  G(t)  in  Eq.  (5-38)  can  be  shown  to  be: 

G(t)  =  ST(t)0(t)S(t)  +  S^tjp-1  S ( t  )  (5-48) 

f  ss  f 

which  can  be  easily  verified  by  direct  differentiation. 

5.3.5  Calculation  of  the  Optimal  Control 

The  control  is  computed  by  writing  Eq.  (5-39)  in  the  form 
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u(t) 


=  -C(t)  x(  t)  -  D(t)t|> 

=  -R-1BT[p  x(t)  +  ©“VjxU)  +  S(t)Xl  (t)  ]  (5-49) 

where  x-|(t)  is  obtained  by  solving  the  following  linear  equation: 

G(t)x^  (t)  =  tj;  -  S(t)x-(t) 


5.3;6  Example  Maneuvers 

Ihe  vehicle  configuration  and  system  parameters  of  Section  5.2.5 
are  assumed  for  the  example  maneuvers  of  this  section.  Furthermore,  full 
state  feedback  is  assumed  for  the  results  of  this  section  and  the  output 
vector  is  assumed  to  be  defined  by  Eq.  (5-31). 

Referring  to  Tables  5-1b  and  5-3,  the  graphical  summaries  of  the 
states  and  controls  are  discussed  qualitatively  in  what  follows. 

Case  1  (Fig.  5-7)  presents  a  terminal  controller  example  maneuver 
in  which  two  flexible  modes  are  controlled  in  addition  to  the  rigid  bod/ 
rotation.  There  are  six  states  in  this  problem,  all  of  which  are 
specified  at  the  final  time.  As  a  result,  P,  S,  and  G  in  Eqs.  (5-42), 
(5-45) ,  and  (5-48)  are  6x6  matrices.  The  time  histories  for  this  case 
are  virtually  identical  to  the  time  histories  obtained  from  the 
corresponding  open-loop  controller,  using  the  same  weight  matrices  in  the 
performance  index.  The  observed  jump  discontinuities  in  the  control 
torques  at  the  initial  and  final  times  are  characteristic  of  this  type  of 
control  design. 

Case  2  (Fig.  5-8)  presents  a  spin-up  maneuver  with  two  flexible 
modes  controlled,  in  which  the  final  angle  is  free  to  be  determined  by 
the  controller.  As  a  result,  the  P,  S,  and  G  matrices  have  the  following 
dimensions:  P(6x6),  S(6x5),  and  G(5x5).  In  particular,  the  final  angle 
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selected  is  exactly  the  angle  selected  by  the  analogous  open-loop  free 
final  angle  transversality  conditions  [Ref.  5-1]. 

Case  3  (Fig.  5-9)  presents  the  same  maneuver  as  in  Case  1,  except 
that  the  performance  index  now  includes  penalties  on  the  first  and  second 
time  derivatives  of  the  control  torques  (Ref.  5-2] .  Including  these 
penalties  allows  the  control  torques  and  torque-rates  to  be  specified  at 
the  initial  and  final  times,  thereby  eliminating  the  terminal  jump 
discontinuities  in  the  control  profiles.  For  this  particular  control 
problem,  the  states  are  augmented  to  include  the  control  and  first  time- 
derivative  of  the  control,  while  the  second  time-derivative  of  the 
control  becomes  the  commanded  input  to  the  augmented  system.  On 
comparing  with  the  results  of  Case  1.  one  finds  a  much  smoother  modal 
amplitude  and  control  torque  time  history,  although  the  peak  torques  and 
tip  deflection  are  somewhat  higher  than  those  in  Case  1.  The  dimensions 
of  the  P,  S,  and  G  matrices  are  all  10x10  for  Case  3. 

Case  4  (Fig.  5-10)  presents  the  same  maneuver  as  in  Case  2,  with 
the  additional  penalty  on  the  control-rates.  As  in  Case  2,  all  the 
states,  controls,  and  control-rates  are  specified  at  the  final  time  for 
the  terminal  controller,  except  for  the  maneuver  angle,  which  is  not 
specified  at  the  final  time.  As  a  result,  the  p,  s,  and  G  matrices  have 
the  following  dimensions:  P(10x10),  S(10x9),  andG(9x9).  On  comparing 
the  results  with  Case  2,  the  modal  amplitude  and  control  torque  time 
histories  are  found  to  be  much  smoother,  although  the  peak  values  are 
higher.  Moreover,  the  terminal  jump  discontinuities  in  the  control 
torque  shewn  in  Case  2  are  eliminated  in  Case  4,  thus  making  it  less 
susceptible  to  control  spillover. 

Case  5  (Fig.  5-11)  presents  the  results  of  a  maneuver  in  which  the 
weight  matrices  and  final  conditions  are  identical  with  those  of  Case  3. 
The  initial  conditions  are  different  from  those  in  Case  3,  but  the 
time-varying  feedback  gains  are  identical.  As  the  plots  in  Figure  5-10 
show,  the  terminal  controller  has  little  difficulty  in  bringing  the 
system  to  the  required  final  conditions. 
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Case  6  (Fig.  5-12)  presents  the  results  of  a  maneuver  in  which  the 
moment  of  inertias  of  the  hub  and  appendages  ''re  increased  by  10%  and  the 
appendage  mass  per  unit  length  is  decreased  by  10%,  while  the 
time-varying  feedba'ck  gains  are  computed  using  the  nominal  structural 
parameters.  Although  the  structural  parameters  are  changed  by  only  10%, 
the  overall  effect  on  the  entire  structure  is  quite  large.  Specifically, 
the  first  and  second  mode  eigenvalues,  u>2,  are  increased  by  16%  and 
22%,  respectively,  while  the  ratio  of  the  hub  moment  of  inertia  to 
appendage  moment  of  inertia  is  increased  by  22%.  However,  the  plots  show 
that  with  the  terminal  controller,  most  of  the  state  variables  reach 
their  prescribed  terminal  boundary  conditions. 

The  results  of  in  Cases  1  through  5  have  been  compared  with  the 
corresponding  open-loop  maneuvers,  and  have  been  found  to  be  identical 
within  plotting  accuracy,  as  expected.  Furthermore,  the  closed  form 
solutions  for  P(t),  S(t),  and  G(t)  have  also  been  compared  with  the 
backward  integrated  solutions  for  verification. 


» 

5.4  An  Analytic  Solution  for  the  State  Trajectories  of  a  Feedback 

Control  System 

As  part  of  the  normal  process  of  a  control  system  design,  the 
analyst  typically  is  interested  in  determining  the  state  trajectories  for 
the  controlled  system.  In  practice,  this  process  is  straightforward, 
since  the  feedback  form  of  the  control  can  be  introduced  in  the  equation 
of  motion  and  numerically  integrated.  Nevertheless,  this  process  can  be 
computationally  intensive  if  either  time-varying  control  gains  are  used 
or  small  integration  step  sir.es  are  required  by  the  presence  of  high 
frequency  system  dynamics. 

In  an  effort  to  overcome  the  computational  difficulties  listed 
above,  a  change  of  variables  is  presented  in  this  subsection  for  the 
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standard  closed-loop  system  dynamics  equation,  which  permits  a  closed 
form  expression  to  be  obtained  for  the  state  trajectories. 


5.4.1  Optimal  Control  Problem 

The  fixed  time  linear  optimal  control  problem  is  formulated  by 
finding  the  control  inputs  u(t)  to  minimize 

t£ 

J  =  '  ~  x‘F  SFxf  +  ~  f  [x  F  QFx  +  u  Ru)  dt  (5-50) 

fco 

for  the  system 

x  =  Ax  +  Bu  ,  given  x(tQ)  (5-51) 

y  =  Fx  (5-52) 

i 

where  x,  u,  A,  B,  F,  Q,  R,  and  S  are  defined  following  iig.  (5-3). 

As  shown  in  Reference  5-6,  the  optimal  control  'a  given  by 

u  ( t )  =<  -R~1BTP(t)x(t)  (5-53) 

where  P  is  the  solution  for  the  differential  matrix  Riccaci  equation 

P  »  -AP  -  PAT  +  PBR'1BTP  -  Q  ;  P(tf)  =  S  (5-54) 

Upon  introducing  Eq.  (5-53)  into  Eq.  (5-51),  the  standard  closed- 
loop  system  dynamics  equation  follows  as 
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x(t)  =  [a  -  BR-1BTp(t)  ]x(t)  ;  x  =  x(t) 


(5-35) 


(It  should  be  observed  that  the  equation  above  is  the  adjoint  to  the 
homogeneous  part  of  Eq.  (5-11).) 

To  obtain  the  solution  for  Eq.  (5-55),  the  following  closed  form 
solution  for  P(t)  (Refs.  5-1,  5-3]  is  introduced  in  Eq.  (5-55): 


P(t)  =  P  +  0  (t) 

ss 


(5-56) 


where  Pgs  is  is  defined  by  Eq.  (5-43).  The  variable  0(t)  is  given  by 


A(t-t  )  AT(t~t  ) 

0(t)  a  0  +  e  |(F  SF-  P  )  -  0  ]e 

SS  t  SS  SS  J 


(5-57) 


-  -IT  ( • ) 

where  A  =  A  -  BR  B  Psg»  e  is  the  exponential  matrix,  and  the  solution 
for  0SS  is  defined  by  Eq.  (5-16). 


Substituting  Eq.  (5-56)  into  Eq.  (5-55),  yields  the  modified  form 
of  the  closed-loop  system  dynamics  equation 


x(t)  =  [a  -  BR  Ve  1  { t)  ]x( t)  /  xq  =  x(tQ)  (5-58) 


where  it  follows  that  the  equation  above  is  nonautonomous. 


5.4.2  Change  of  Variables  for  x(t) 


To  simplify  Eq.  (5-58)  the  following  coordinate  transformation  for 
the  dependent  variable  x(t)  is  introduced: 


x(t) 


0(t)r(t) 


(5-59) 


where  0(t)  is  given  by  Eq.  (5-57)  and  r(t)  is  a  vector  function  which  is 
to  be  determined. 

Upon  differentiating  Eq.  (5-59)  we  find 


x  =  0r  +  0r 


(5-eo) 


or 


-T  -1  T\  • 

=  (A0  +  0A  -  BR  B  Jr  +  0r 


(5-61) 


.  .  _  _T  - 1 

where  0  in  Eq.  (5-60)  has  been  replaced  by:  0  =  A0  +  0A  -  BR  B. 

The  differential  equation  for  r  is  obtained  by  introducing  Eqs. 
(5-59)  and  (5-61)  into  Eq.  (5-56),  leading  to 


0(r  +  ATr)  *  0  (5-62) 

from  which  it  follows  that  the  linear  constant  coefficient  vector 
differential  equation  for  r  is  given  by 


•  rT 

r  =  -A  r 


r 

o 


0~V  ) x 

o  o 


(5-63) 


The  solution  for  r  follows  as 


-T 

-A  (t-t  ) 

r(t)  =  e  °  r  (5-64) 

o 

Substituting  Eq.  (5-64)  into  Eq.  (5-59)  produces  the  desired 
solution  for  the  state  trajectories  as 


x(t) 


(5-65) 


-T 

-A^t-t  ) 

=  9(t)e  °  r 


5.4.3  Recursion  Relationship  for  Evaluating  the  State  at  Discrete  Times 

If  the  solution  for  x(t)  is  required  at  the  discrete  times 
31  fco  +  kAt  (k  =>  1 ,  •  •  «,N)  for  At  =  (tf  -  t0)/N,  then  Eq.  (5-65) 
can  be  written  as 


X<V  •  V-  +  \  •  k  ■ 


(5-66) 


where 


-T 

At 


A  =»  0  e 

O  S3 


-At.  n,  ,  AT(t  -t. 

f r , .JT_ _  _  .-1  .  .  1  of 

e  (F  S?  -  P  )  -  0  e  i 

L  ss  ssJ 


Ak  a  Ak-ie  ^  »  b^  =  e  ^^b^  ^  i  k=;,*»»,N 


5.4.4  Conclusions 

Closed  form  solutions  for  the  feedback  gains  required  by  an 
optimal  linear  tracking  controller  and  an  optimal  terminal  controller 
have  been  developed.  Results  of  example  maneuvers  have  been  shown  which 
demonstrate  the  efficiency  and  validity  of  the  formulations  described. 

The  use  of  control-rate  penalties  in  each  case  have  been  shown  to  improve 
the  overall  system  response. 
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In  addition,  a  straightforward  algorithm  has  been  presented  for 
generating  the  state  trajectories  for  a  feedback  control  system.  The 
algorithm  is  computationally  efficient  in  that  no  numerical  integration 
is  required  and  simple  recursion  relationships  generate  the  desired 
solution  at  discrete  times.  Furthermore,  this  algorithm  has  significant 
potential  if  used  in  conjunction  with  algorithms  which  attempt  to  enhance 
system  robustness,  by  iteratively  refining  the  weighting  matrices 
appearing  in  Eq.  (5-50). 
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No.  of 
Modes 

t  -t 
f  o 

(sec) 

*(1) 

8 

0 

(rad) 

• 

e 

o 

(i*ad/s ) 

9f 

(rad) 

• 

9f 

(rad/s) 

CPU(2) 

(sec) 

2 

15 

0 

0 

0 

0.4 

0 

27 

2 

15 

2 

0 

0 

0.4 

0 

•59 

2 

15 

2 

0.1  5 

-0.02 

0.4 

0 

58 

2 

1.5 

2 

0 

0 

0.4 

0 

58 

Case 

No. 

No.  of 
Modes 

t,-t 
f  o 

(sec) 

k("  9o 
(rad) 

• 

9 

o 

(rad/s  ] 

9f 

(rad) 

»* 

(rad/s ) 

CPU(2) 

(sec) 

1 

2 

15 

0 

0 

0 

0.4 

0 

34 

2 

2 

15 

0 

0 

0 

(free) 

0.04 

32 

3 

2 

15 

2 

0 

0 

0.4 

0 

83 

4 

2 

1  5 

2 

0 

0 

(free) 

0.04 

77 

5 

2 

15 

2 

0.15 

-0.02 

0.4 

0 

83 

6 

2 

15 

2 

0 

0 

0.4 

0 

83 

Table  5-2.  Weighting  matrices  for  example  maneuvers: 
linear  tracking  problem. 


Case  No. 

Open-Loop  Weighting  Matrices^ 1 • 

1 

Q*  =  diag  [l(-5)  1(-3)  1<-3)  1  (-3)  K-3)  1  (-3)  ] 

R*  «  diag  [  1 (-3)  1(0)  ] 

1 

Q*  =  ®  W  ®  W<4 

EE  oo  11 

W  »  diag  [l(-5)  1  (-3)  1  (-3)  1  (-3)  K-3)  1(~3)  ] 

Wqo  =  diag  [l  (-9)  K-9)  ] 

W  -  diag  [l (-9)  1(-9)  ] 

R*  =  diag  (l  (-3)  1(1)] 

Closed-Loop  Weighting  Matrices 

1 

Q  =  diag  [l (-3)  1(-2>  1(-3)  1 (-5)  ] 

R  =  diag  [l (-3)  1 (0)  ] 

n 

Q  -  diag  (l(-3)  1(-2)  1.(-3)  1(-5)  1(-9)  K-9)  1(-9)  1  (-9)  ] 
R  =  diag  [l(-3)  1(1)  ] 

Notes: 


1.  b(a)  denotes  b  x  1 0a 

2.  ®  denotes  direct  sum 


b(a)  denotes  b  x.  1 0a 


TIME  III 


15 


0.40 


Figure  5-11. 


Case  5,  spin-to-rest  maneuver,  2  modes,  5  controls, 
control-rate  penalty  with  k  «*  2. 
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